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DYNAMICS OF THE M87 GLOBULAR CLUSTER SYSTEM 



ABSTRACT 

Wc present a dynamical analysis of the globular cluster system associated with M87 (= NGC 4486), 
the cD galaxy near the dynamical center of the Virgo cluster. The analysis utilizes a new spectroscopic 
and photometric database which is described in a companion paper (Hanes et al. 2001). Using a sample 
of 278 globular clusters with measured radial velocities and metallicitics, and new surface density profiles 
based on wide-field Washington photometry, we study the dynamics of the M87 globular cluster system 
both globally — for the entire cluster sample — and separately — for the metal-rich and metal-poor 
globular cluster samples. This constitutes the largest sample of radial velocities for pure Population II 
tracers yet assembled for any galaxy. Our principal findings are summarized as follows: 

(1) Surface density profiles constructed from our Washington photometry reveal the metal-poor clus- 
ter system to be more spatially extended than its metal-rich counterpart, consistent with earlier 
findings based on HST imaging in the central regions of the galaxy. Beyond a radius of i? ~ 1.5Re 
(10 kpc), the metal-poor component dominates the total globular cluster system. 

(2) When considered in their entirety, each of the combined, metal-poor and metal-rich globular cluster 
samples (278, 161 and 117 clusters, respectively) appear to rotate, with similar rotation amplitudes, 
about axes whose position angles are indistinguishable from that of the photometric minor axis, 
Oo = 65°. 

(3) The one-dimensional rotation curve (i.e., binned in circular annuli) for the metal-rich cluster system 
has a roughly constant mean amplitude of flR = IGO^gg km s^^. The metal-rich clustcirs appear 
to be rotating, at all radii, about the photometric minor axis of the galaxy. However, a smoothed, 
two-dimensional map of the line-of-sight velocity residuals suggests that the rotation field for the 

metal-rich clusters is non-cylindrical in nature. Instead, it exhibits a "double-lobed" pattern, with 
maxima at i? 3.5-4i?e (25-30 kpc) along the approximate photometric major axis of the galaxy. 

(4) The one-dimensional rotation curve of the metal-poor cluster system has mean amplitude of flR = 
1721^03 km s~^. The two-dimensional map of the rotation field for the metal-poor clusters shows 

some evidence for solid-body rotation or, alternatively, for a "shear" in the line-of-sight velocity. 
This shear is similar in size and orientation to that observed for Virgo galaxies within two degrees of 
M87, and is consistent with a scenario, previously suggested on the basis of dwarf galaxy kinematics 
and x-ray imaging, in which material is gradually infalling onto M87 along the so-called "principal 
axis" of the Virgo cluster. 

(5) Beyond a radius of i? 2Re (15 kpc), the approximate onset of the galaxy's cD envelope, the metal- 
poor globular cluster system rotates about the photometric minor axis similar to its metal-rich 
counterpart. Inside this radius, however, the metal-poor clusters appear to rotate around the 
photometric major axis. 

(6) The complete sample of 278 globular clusters has an almost perfectly isotropic velocity ellipsoid, 
with = 1 - CT^/cr^ ~ 0. 

(7) When considered separately, the metal-poor cluster system shows a modest but significant tan- 
gential bias of Pc\ — —0.4, while the velocity ellipsoid of the metal-rich cluster system is radially 
biased, with (id — -1-0.4. 

Taken together, these results demonstrate that the dual nature of the M87 globular cluster system 
— first identified on the basis of its bimodal metallicity distribution - - also extends to its dynamical 
properties. We discuss the implications of these findings for the various formation scenarios proposed 
for giant elliptical galaxies and their globular cluster systems. 

Subject headings: galaxies: halos galaxies: clusters galaxies: individual (M87) — galaxies: 
kinematics and dynamics - galaxies: star clusters 

1. INTRODUCTION 

As surviving relics from the epoch of galaxy formation, globular clusters (GCs) provide important insights into the 
processes which shaped their present day host galaxies. The obscirvcid properties of GCh in giant elliptic;al galaxies such 
as their total numbers, chemical abundances and spatial distributions — have been used alternatively to argue for a wide 
diversity of formation processes including dissipative (monolithic) collapse, spiral-spiral mergers (Ashman & Zepf 1992), 
multi-modal star formation histories caused by galactic winds (Harris, Harris & McLaughlin 1998), and dissipationless 
hierarchical growth (Cote, Marzke & West 1998; Cote et al. 2000). There is, however, one important piece of observational 
evidence which has been almost entirely absent in such discussions — the dynamical properties of the GCh. Due to the 
faintness of even the most luminous GCs in nearby giant elliptical galaxies, the requisite spectroscopic observations are 
extremely challenging, and it is only recently that radial velocities for significant samples of clusters have been accumulated 
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in a handful of galaxies (e.g., Cohen & Ryzhov 1997; Sharpies et al. 1998; Zepf et al. 2000). 

M87, the incipient cD galaxy near the dynamical center of the Virgo cluster, offers an ideal testing ground for the 
various formation scenarios. It has by far the most populous globular cluster system (CCS) of any nearby galaxy (e.g., 
Nci ~ 13,500; McLaughlin et al. 1994), and a wealth of observational data has been accumulated for both the galaxy 
itself, and the surrounding Virgo cluster. Such material includes numerous analyses of the photometric properties and 
spatial distribution of M87 GCs {e.g., Strom et al. 1981; McLaughlin et al. 1994; Geisler, Lee & Kim 2001), surface 
brightness profiles for the galaxy light (de Vaucouleurs & Nieto 1978; Carter & Dixon 1978), stellar dynamical studies 
of the central regions of M87 (van der Marel 1994), photometric and radial velocity surveys of surrounding Virgo cluster 
galaxies (Binggeli, Sandage & Tammann 1985; Binggeli, Popescu & Tammann 1993; Girardi et al. 1996), and X-ray 
observations of the hot intracluster gas which fills the Virgo gravitational potential well (Nulsen & Bohringer 1995). 

The dynamics of the M87 GCS have been addressed in several previous studies, albeit in varying levels of detail. Early 
efforts (Huchra & Brodie 1987; Mould, Oke & Nemec 1987; Mould et al. 1990) were hampered by limited sample sizes: 
for example, the compilation of Mould et al. (1990) consisted of radial velocities for only 43 clusters. Although these data 
were sufhcient to provide a provisional mass estimate for M87, and to hint at the presence of rotation in the GCS, it was 
shown by Merritt & Tremblay (1993) that this sample provides only weak constraints on the distribution of dark matter 
in the M87 halo. Moreover, Merritt & Tremblay (1993) demonstrated that a measurement of the velocity anisotropy of 
GCS would require radial velocities for a minimum sample of several hundred clusters. 

In a major observational effort, Cohen & Ryzhov (1997) used the Keck telescope to measure radial velocities for a 
total of 205 GCs surrounding M87. Using these data, Cohen & Ryzhov (1997) found clear evidence for a rising velocity 
dispersion profile, and derived a total mass of Al = (3.8 ± 1.2) x lO^^M© within a volume of radius of r ~ 33 kpc centered 
on M87. An independent analysis of these velocities by Kissler-Patig & Gebhardt (1998) suggested that the M87 GCS 
is rotating rapidly at large radii: i.e., flR ~ 300 km s~^ beyond a projected radius of R — 3i?e (20 kpc)]^ This rapid 
rotation was subsequently confirmed by Cohen (2000) who augmented her original radial velocity sample to include 17 
additional clusters having R > 3i?e. 

While the amount of angular momentum implied by these results is large, there is emerging evidence from radial velocity 
surveys of planetary nebulae that giant elliptical galaxies may, as a class, rotate rapidly at large radii (e.g., Hui et al. 1995; 
Arnaboldi et al. 1996; 1998). An enormous amount of theoretical effort has been devoted to understanding the acquisition 
of angular momentum during galaxy formation: e.g., via tidal torques in monolithic collapse scenarios (Peebles 1969), 
or through the conversion of orbital and spin angular momentum of the progenitor galaxies in major (Heyl, Hernquist 
& Spergel 1996) and multiple (Weil & Hernquist 1996) mergers. The angular momentum content of elliptical galaxies 
may thus constitute an important clue to their formation; ultimately, the rotation properties of the chemically distinct 
GC populations in these galaxies may provide a powerful means of discriminating between the various formation models 
described above. 

In this paper, we present a dynamical analysis of the M87 GCS using an improved and expanded sample of radial 
velocities. In addition, new Washington CTi photometry is used to derive metallicities for all program objects, and to 
investigate correlations between metallicity and dynamics by carrying out separate analyses for both the metal-rich and 
metal-poor GC populations. 

2. THE DATABASE 

The database used in the analysis presented below is a subset of that which has been described in detail in a companion 
paper (Hanes et al. 2001; hereafter Paper I), so only a brief summary is given here. The reader is referred to Paper I for 
a complete description of the sample properties, observational material and data reductions. 

2.1. Photometry and Radial Velocities 

One of our goals is to investigate the relationship between GC dynamics and metallicity, so we require homogeneous 
metallicity estimates for our sample of clusters. Such information is most easily obtained from broadband photometry. 
Unfortunately, although there have been numerous photometric studies of the M87 GCs, the large areal coverage of our 
radial velocity survey means that no single existing survey includes all of our program objects. For instance, the most 
extensive publsihed photometric survey of the M87 GCS remains that of Strom et al. (1981), who presented photographic 
UBR photometry for 1728 candidate GCs surrounding M87. Yet our sample includes many objects not included in the 
Strom et al. (1981) catalog: e.g., objects that are at large radii, beyond the limits of their survey, or close to the galaxy 
center where crowding and background contamination is severe. Therefore, as described in Paper I, we have obtained 
deep, wide-field CCD images centered on M87 using the CTIO and KPNO 4m telescope. Imaging was carried out in 
Washington C and Ti filters, since the long color baseline of the C— Ti index provides exceptional metallicity sensitivity 
(see, e.g., Geisler, Lee & Kim 1996). The photometric observations, calibrations, and reductions have been described 
briefly in Paper I, while a more comprehensive description of the photometric properties of the M87 GCS based on these 
data will be presented in a forthcoming paper (Geisler, Lee & Kim 2001). 

The sample of radial velocities we begin with here consists of data for 334 candidate and previously confirmed GCs. 
This number includes 145 measurements from the Canada- France- Hawaii telescope (CFHT), 87 of which are for targets 
not observed in any of the previous surveys of Huchra & Brodie (1987), Mould, Oke & Nemec (1987), Mould et al. (1990), 

^^According to de Vaucouleurs & Nieto (1978), iJe = 96" for M87. Throughout this paper, we adopt an M87 distance of D = 15 Mpc (Pierce 
et al. 1994; Ferrarese et al. 1996), so that 1' = 4.364 kpc. 
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Fig. 1. — Upper Left Panel: Distribution of (C— Ti) colors for 334 globular cluster candidates with measured radial velocities (solid 
histogram). For comparison, the color distribution of all 2744 globular cluster candidates from the photometric study of Geisler, Lee & Kim 
(2001) is shown as the dashed histogram. The dashed vertical lines indicate our adopted selection criteria on color, 0.80 < (C— Ti) < 2.35 
mag, while the dotted vertical line shows the color, (C— Ti) = 1.42, used to isolate the blue and red cluster samples. Upper Right: Radial 
velocity histogram for the 334 candidate globular clusters with accurate (C— Ti) colors. The vertical lines show our adopted selection criteria 
on velocity: 200 < v < 2550 km s~^. Lower Panel: Color vs. radial velocity for our initial sample of 334 objects. The dotted region shows 
the joint selection criteria on color and radial velocity. This box defines a final sample of 278 bona fide globular clusters. 

Cohen & Ryzhov (1997), or Cohen (2000). The total is down slightly from the full database of 352 objects in Paper I, 
as we have discarded out of hand 12 targets for which we were unable to obtain unambiguous Washington photometry 
(ID nos. 881, 5053, 5055, 5058, 5064, 5065, 5066, 5067, 5071, 9002, 9051, and 9052 from Paper 1); five which were already 
identified by Cohen & Ryzhov (1997; CR97) as Galactic M stars (ID nos. 38, 87, 782, 1551, 8001); and one object which 
turns out to be visibly non-stellar (ID no. 207). 

All measurements have been transformed onto common astrometric and radial velocity scales. We adopt the radial 
velocity and uncertainty from our own analysis for all targets observed at the CFHT and transform all other data to this 
system as described in Paper I. Briefly, cluster candidates that were not observed at the CFHT, but were observed at 
Keck by Cohen & Ryzhov (1997), are transformed onto our MOS velocity scale according to: 

Up(MOS) = 0.96wp(CR97) + 91 

and assigned an uncertainty of ±100 km s~^, as recommended by these authors.0 Objects measured neither by us nor 
by Cohen & Ryzhov (1997), but by Cohen (2000), are transformed as: 

t;p(MOS) = 0.98wp(C2000) + 120 

and given an uncertainty of ±50 km s~^. Finally, to targets with velocities recorded only by Mould et al. (1987) or Mould 
et al. (1990), we apply the correction: 

Wp(MOS) = 0.76Dp(Mould) + 280 

and adopt an uncertainty of ±200 km s~^. 

2.2. Sample Selection 

Before proceeding with the dynamical analysis, it is important to first identify a sample of bona fide M87 GCs. Culling 
Galactic field stars from the sample on the basis of radial velocity is not trivial since the relatively low systemic velocity 
of M87 (1277 ± 2 km s^^; van der Marel 1994) and the high velocity dispersion of the GCs (i.e., (ap) ~ 400 km s~^ see 
§3) make the correct choice of the low velocity cutoff problematical. We adopted a selection criterion of: 

200 < Dp < 2550 km s"\ (1) 

^^Here and throughout, quantities referred to the line of sight are denoted by a subscript p. 
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Fig. 2. — Spatial distribution of the 278 confirmed globular clusters around M87 (i.e., all objects falling within the boxed region of the 
bottom panel of Figure hh. The center of the galaxy is at (0,0), the inner circle marks the effective radius of the stellar mass of M87 
{Re = 96" ~ 7 kpc), and the outer circle has a radius of 414" ~ 30kpc ~ 4.3-Re. The diagonal line shows the photometric major axis of the 
galaxy. 



This range is roughly symmetric about the mean velocity, (vp) = 1350 km s~^, of the M87 GCS (see §3) and is expected 
to exclude the vast majority of foreground stars and background galaxies. We also imposed a selection on C— Ti color, 
discarding all objects which did not have colors within the range: 

0.8 < C - Ti < 2.35 mag. (2) 

The Galactic foreground extinction in the direction of M87, according to the DIRBE maps of Schlegel, Finkbeiner & 
Davis (1998), is E(B-V) = 0.022 mag, which corresponds to E(C-Ti) = 0.045 mag (Seeker et al. 1995). The metalhcity 
calibration of Geisler & Forte (1990) suggests that this color selection will confine our sample to clusters having metallicities 
in range —2.6 < [Fe/H] < 1.0 dex, an interval which should encompass virtually all GCs in M87 (e.g., Cohen, Blakeslee & 
Ryzhov 1998). These selection criteria are illustrated in the upper panels of Figure |. The filled histo grams in these two 
panels show the color and velocity distributions for all 334 candidate GCs with measured radial velocities and (C— Ti) 
colors. The dashed vertical lines in each panel show the adopted upper and lower selection limits. The dotted vertical 
line in the upper left panel shows our adopted dividing point between the metal-rich and metal-poor cluster populations: 
C— Ti — 1.42. This color, which corresponds to a metalhcity of [Fe/H] = —1.16 dex, was chosen based on the observed dip 
in the (C— Ti) distribution of the 2744 globular cluster candidates from the photometric catalog of Geisler et al. (2001) 
which have Ti < 22 (i.e., the dashed histogram in this panel). 

The lower panel of Figure Q is a color-velocity diagram for these same objects. The dotted region shows the joint 
selection criteria on color and velocity, which together identify a total of 278 bona fide GCs. Of these, 161 are classified 
as "blue" (C— Ti < 1.42), or metal-poor, and 117 are considered "red" (C— Ti > 1.42), or metal-rich, for purposes of the 
analysis and modelling in §§3 and 4 below. Note that three objects (ID nos. 161, 508, and 1602) have radial velocities 
which are appropriate for M87 GCs, and yet have C— Ti colors which are far too red for true clusters (i.e., the inferred 
metallicities are typically several hundred times above solar). A fourth candidate (ID no. 5002) has a radial velocity near 
the mean of the GC sample, but an anomalously blue color. Although we have no reason to suspect either the colors or the 
velocities of these objects, their nature is unclear and additional observations are clearly warranted, first and foremost to 
confirm our results (particularly for Strom 508, an unresolved object which has a high radial velocity of Vp = 5817 ± 120 
km s~^ but a very red color of C— Ti ~ 3.0). In any event, all four are omitted from the following analysis. 

Our innermost GC is projected to a distance of roughly 36" from the center of M87 (or about 2.6 kpc for our assumed 
distance of 15 Mpc to the galaxy), while the outermost member of our sample lies at i? ~ 631" ~ 46 kpc. Figure || shows 
the distribution on the sky of all 278 bona fide GCs (i.e., those objects falling inside the dotted region in the lower panel 
of Figure |l|). The solid points on this plot denote 109 clusters which were observed by us at the CFHT; of these, 54 
are entirely new measurements and the rest are repeat observations of GCs identified in one of the earlier studies. The 
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open symbols refer to data taken from these other surveys, as indicated in the legend. The straight line running through 
the figure represents the photometric major axis of the M87 halo light (oriented at 155° East of North; Carter & Dixon 
1978), the inner circle marks the effective radius of the galaxy light {Re = 96" ~ 7 kpc [de Vaucouleurs & Nieto 1978]), 
and the larger, dashed circle has a radius of 414" ~ 30 kpc ~ 4.3i?e. Outside this circle, the azimuthal coverage of the 
composite radial-velocity database is far from complete, a situation which can significantly compromise any inferences on 
the kinematics of the CCS at such large galactocentric distances (see below). 

3. KINEMATICS OF THE CLUSTER SYSTEM 

In previous analyses of the kinematics of the M87 CCS, investigators have usually proceeded by fitting sine curves to 
their line-of-sight velocities, Up, as a function of projected azimuth, O, in order to test for the presence of rotation and 
to measure its amplitude (see Cohen & Ryzhov 1997; Kissler-Patig & Gebhardt 1998; Cohen 2000) We shall do the same 
here, but we pause first to clarify precisely what this assumes about the intrinsic velocity field of the GCS. 

3.1. Mathematical Context 

Consider a spherical GCS rotating about some axis, with an angular velocity oj that may be any axisymmetric function 
of position inside the galaxy.[^ Let r, 6, and (j) denote the usual spherical coordinate system within the GCS, such that 
^ = naturally identifies its rotation axis and its volume density n is a function only of the galactocentric distance r. 
Also define a projected, circular coordinate system on the plane of the sky by the radial coordinate R and the azimuth 
0. The projected (surface) density of the GCS, N, is, by the spherical hypothesis, a function of R alone: 

f°° rdr 

N{R) = 2 n{r)^==. (3) 

However, its projected angular velocity, 17, which is just the average of lo along the line of sight through the GCS, may be 
a function of R and O both. Assuming that the cluster system is viewed so that its rotation axis lies exactly in the plane 
of the sky and coincides with O = 0, we have 

where cos^ — (_Rcos8)/r by simple geometry (see also Fillmore 1986). In addition to this, the angular velocities uj and 
n are related to the linear velocities {v^) — the intrinsic streaming, or rotation, velocity of the GCSF^ — and (up) — 
the averaged component of velocity along the line of sight — in the usual way: lo = {vij,) / {r snvO) and 12 = (z)p)/i?sinO. 
Thus, 

NiR) A "(^Vsin^ V^^31?2 ' 

where (v^) may still, in general, be any function of r and 9. 

Equation (^ makes it clear that the apparent angular velocity will be a function of R alone, fl ~ ^{R), if ti^ is a function 
only of r, i.e., if the intrinsic angular velocity of the GCS is constant on spherical surfaces. This is just what is required, 
as equation (ra) shows explicitly, to have (vp) cx sinO, and thus to justify fitting a sinusoid to our radial velocity data as 
a function of projected azimuth. 

It is worth emphasizing this point, even though it is not new (see, e.g., Fillmore 1986): the projected angular velocity 
of a spherical GCS is constant on circles in the plane of the sky — so that its linear velocity (vp) varies sinusoidally 
with projected azimuth O — if the intrinsic angular velocity is constant on spheres, such that (v^) cx sinO. It happens 
that the variation of (vp) with O in our sample of M87 GCs is adequately described, at all galactocentric radii, by a 
simple sine curve (see below), and we interpret this as an indication that the angular velocity field of the GCS is stratified 
approximately on spherical surfaces: (v^) / {r sm9) = uj{r). This interpretation is somewhat different from that adopted by 
Cohen & Ryzhov (1997), and apparently by Kissler-Patig & Gebhardt (1998), who fit sines to their Vp vs. O data and then 
posit that the intrinsic {vt/,) field of the GCS is stratified on cylinders. The only way that both of these assumptions can 
be simultaneously true — that is, (vp) oc sinO and (v^) depending on r and 9 only through the combination (r sin0) — is 
in the special case of solid-body rotation, uj = = constant. This would, of course, require that flR oc R. Unfortunately, 
the data neither rule out this possibility nor present a compelling case for it (see, e.g.. Figure |8| below) . We shall return 
to the issue of possible solid-body rotation below, and in §5.7 for the specific case of the metal-poor GCS (see Figure . 

In summary, we assume that: (1) the GCS of M87 can be approximated as spherical [n = n(r) and N = N{R)]; (2)its 
angular velocity field is constant on spheres [uj = uj{r) and fl = fl{R)]; and (3) its rotation axis lies exactly in the plane of 
the sky and coincides with G = 0. We are then able to investigate the rotation of the GCS by fitting our observed radial 
velocities with the function 

Vp{Q) = Vsys + (nR) sin(e - Go) , (6) 

where Oq (measured, as O is, in degrees East of North), locates the GCS rotation axis, and where the rotation amplitude 
flR can in principle be any function of projected galactocentric radius R. 

^^While it is clearly imperfect (see, e.g., McLaughlin, Harris, & Hanes 1994), the assumption of spherical symmetry in the M87 GCS is not 
grossly in error, and we invoke it as a simplifying approximation throughout this paper. 
^"The (/> direction wraps around the rotation axis by definition. 
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3.2. Global Kinematic Properties 

Figure ^ shows the dependence of radial velocity on azimuth for our full sample of 278 GCs in M87, for our subsample 
of 161 blue clusters (defined by 0.80 < C-Ti < 1.42), and for our subsample of 117 red GCs (1.42 < C-Ti < 2.35). 
Overlaid on the data in each case is the best-fit sine. For the full sample, we find a rotation amplitude of fii? = 169lgy km 
s""'^, where the quoted errors refer to 90% confidence intervals. This rotation amplitude is somewhat larger than the value 
of fii? ~ 100 km s~^ estimated by Cohen & Ryzhov (1997), although still consistent given the rather large uncertainties. 
The best-fit position angle of the rotation axis for the full sample is Oq — 66° ± 35° E of N. The photometric minor axis of 
the underlying galaxy hght is similarly oriented to Ophot = 65° (Carter & Dixon 1978), a correspondence that is consistent 
with our simplifying assumption that we are viewing the GCS essentially edge-on. We separately find flR ~ 172^^Jg km 
s~^ and Oo = 59° ± 52° for the entire sample of metal-poor clusters, and fii? ~ 160l;gg° km s~^ and 8o = 76° ± 45° for 
the metal-rich sample. The dispersion of velocities about the best-fitting sine is cTp^r ^ 384:1:11 km s-^ for the full sample, 
Cp,!- = 397l4g km s~^ for the metal-poor clusters, and CTp^r — 364^52 for the metal-rich clusters. 

We conclude, therefore, that rotation is present in the M87 GCS and is significant at a high confidence level, for the blue 
and red subsystems alike. Moreover, the metal-poor and metal-rich clusters in our sample have: (1) essentially identical 
rotation amplitudes; (2) rotation axes that coincide with one another and with the photometric minor axis of the galaxy; 
and (3) similar fractions of their total kinetic energy stored in rotation [{^lR/ap^Y)^ ~ 0.2 ±0.1 in each case]. Thus, within 
the observational uncertainties, there is no obvious difference between the global kinematics of the two subsystems. We 
caution, however, that this is just a thumbnail sketch of the kinematics, averaged over the full range of galactocentric 
radii spanned by the GCs in our radial-velocity sample. As was stressed above, the quantity QR can, in principle, be a 
function of projected radius — as might be the rotation axis Qq and the dispersion of the line-of-sight velocities. While 
Figure ^ is useful as a rough summary of the situation, it is necessarily limited. 

Thus, we sort our GCs into a number of distinct circular annuli and fit sine curves separately to Vp vs. O in each 
annulus to investigate possible radial variations in the kinematics of the total, blue, and red cluster systems. Table ^ 
lists the results of this exercise and includes full details of the global kinematics discussed just above. The first column 
of the Table defines the annulus in question (with 1" = 0.073 kpc for a distance of 15 Mpc to M87), the second and third 
columns gives the median radius (in arcseconds and kpc) of the clusters in that radial bin, and the fourth records the 
number of objects in it. The fifth and six columns then give the mean line-of-sight velocity of the clusters (ignoring any 
rotation that might be present, and estimated using the robust, biweight measure of "location" discussed by Beers, Flynn, 
& Gebhardt 1990) and the dispersion of their velocities about this mean value (i.e., the biweight "scale" of Beers et al. 
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Fig. 3. — Radial velocity vs. azimuth for the full sample of 278 confirmed globular clusters, and separately for the 161 blue clusters with 
0.80 < (C-Ti) < 1.42 and the Hftred globular clusters with 1.42 < (C-Ti) < 2.35. The best-fit sine curve, Vp = v^ys + (Oi?) sin(e - 60), 
is overlaid in each case; see Table hi for a listing of the parameters. The orientation of the photometric minor axis of M87 is indicated by the 
vertical arrows in each panel at = 65° and 245° East of North. 
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Fig. 4. — Line-of-sight velocity vs. azimuth for M87 globular clusters in six radial bins. From top to bottom, the panels show clusters in the 
range: (1) 35" < R< 125"; (2) 125" < R< 207"; (3) 207" < i? < 297"; (4) 297" < i? < 380"; (5) 380" < i? < 414"; (6) 414" < R< 635". 
Filled circles are metal-poor (blue) globular clusters, and solid lines are their best-fit sine curves with parameters given in Table hi. O 
circles denote the metal-rich (red) globular clusters, and broken hues trace the best-fit sine curves (with parameters also given in Table 
Note that the identical sine curves shown in the bottom two panels have been derived from the combined data: i.e., using all clusters in 
range 380" < H < 635". 



1990). The uncertainties on these values, and on all quantities in Table ll define a 90% confidence interval. They are 
estimated using a numerical bootstrap procedure in which 1000 artificial datasets are constructed by randomly choosing 
N clusters from among the actual data in an annulus. The mean and dispersion are computed for each of these 1000 trial 
datasets, the results are sorted, and the values corresponding to the 5th and 95th percentiles are identified; uncertainties 
are then defined as the offsets between these values and the actual (vp) and cTp computed from the read data. Columns 7, 
8, and 9 of Table ^ give the rotation parameters in each annulus, for each of our GC subsamples. These are determined 
by an error- weighted, nonlinear, least-squares fit of the sinusoid in equation (m, and their 90% confidence intervals follow 
from fitting Vp vs. for each of the 1000 artificial datasets constructed in the oootstrap procedure just described. Finally, 
column 10 of Table |l| gives the (biweight) dispersion of the velocities about the best-fit sine in each annulus for comparison 
with the previous analyses of Cohen & Ryzhov (1997) and Kissler-Patig & Gebhardt (1998). 

Figure directly compares some of the fits from Table |^ against the data. The filled circles in this figure plot Vp 
against ©for our blue GCs, while the open circles correspond to the red clusters. The top panel shows the objects with 
35" < R< 125"; the second panel is for those with 125" < R< 207"; the next panel includes clusters at 207" < R< 297"; 
the fourth is for GCs with 297" < R < 380"; the fifth has 380" < R < 414"; and the sixth, 414" < R < 635". In the 
top four panels, the solid curves trace the corresponding sine fits for the blue clusters, as summarized in Table |^, and 
the broken curves are the best-fit sinusoids for the red clusters. In both of the bottom two panels, the blue and red fits 
shown are those given in Table ^for the annuli combined (i.e., for 380" < R < 635"), as these are much more stable and 
less uncertain. In all cases, the observed fp(0) data appear to be consistent with the simple sinusoidal form that was 
assumed in our fitting. 

Table |l| and Fi gurc H hint at some interesting radial trends in the kinematics of the M87 GCS. First, it appears that 
the projected kinematics of the red and blue subsystems are, at least to first order, rather similar. There are, however, 
suggestions of some differences which we shall return to below and in §5.7 where we discuss the two-dimensional velocity 
field of the GCS. Second, although it is roughly constant over much of our range in projected radius, the velocity dispersion 
of the clusters (either about the sample mean (vp) or about the best-fit sine curve) shows some evidence for an increase 
at the largest R. This is just significant at the 90% level in our sample (and not even then for the metal-rich subsystem), 
but it is actually predicted by the simple dynamical models we present in §4. Third, there is some evidence for a rise 
in the rotation amplitude Oi? at large radii, although this is not formally significant at the 90% level; to this degree of 
accuracy, the red and blue cluster samples both are consistent with a constant flR ^ 170 km s^^ at all radii 2.5 ^ ^ 35 
kpc (see Figure |8| below). On the other hand, the angle of the rotation axis, Oq, does change significantly as a function 
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Table 1 

Kinematics of the M87 Globular Cluster System 



R 

(arcsec) 


(arcsec) 


{R) 
(kpc) 


N 


(km s-i) 


CTp 

(km s-i) 


^sys 

(km s-i) 


eo 

(° E of N) 


VlR 
(km s~^) 


^P>r 

(km s-i) 






Full sample: 


278 clusters with 0.80 < (C - Ti) < 2.35 






oo yjoo 


252 


18.3 


278 


1335t 1? 


40lt f. 


13541 % 


cc+36 
""-34 


iDy_ 97 


■^84+ 
004_ 32 


oo ±ZjO 

125-207 
207-297 
297-380 
380-635 


91 
178 
252 
332 
415 


6.6 
12.9 
18.3 
24.1 
30.2 


55 
57 
55 
56 
55 


1392^ 
1282_ QQ 

1269ti«l 

13111 

1445111° 


3181 % 

■mt i 

409+ % 
3321 ^3 
5431 i 


13791 g 

1312lif4 

12481^29 

13751^^2 
1360llf^ 


02+62 

-19. 

-1+56 
-■^-46 

127l« 
49l« 

00+36 
°^-69 


qn+130 

133 
153_i44 

1 QQ + 138 
100_122 

iqq+113 

OUO_450 


306+ 
3591 ^8^ 

000+ 60 

6m_ 82 
3011 ^9 

4701 


414-635 


402 
481 


29.2 
35.0 


27 
28 


132411°^ 
154711^? 


6051^^6 
460lif4 


1 "51 q+232 
13811^9^6^ 


451^^ 
1051^5 


of;i+436 
o\jJ- 221 

41 1 +305 
^■'■■'■-594 


^14+105 
^-'-^-173 
or9+174 
OOZ_4ig 




Metal-poor 


sample: 161 clusters with 0.80 < (C - Ti) < 1.42 






OO Doo 


260 


18.9 


161 


1339t 


4121 


13751 ^1 


591^3 


1 79+ 51 


007+ 37 


oo L^O 

125-207 

207 297 
297-380 
380-635 


98 
185 
250 
334 
414 


7.1 
13.5 
18.2 
24.3 
30.1 


26 
33 
33 
31 
38 


1355li^6 

1235111^ 

1325ti^4 

1433t}^^ 


3001 

orn+ 96 
OOU_iQ3 

447li«3 
340+ '™ 

r,->-|+ 81 


1277+1^^5 

1318l'9'3° 
144514^4 
1346litl 


77+76 
' '-75 
Q7+58 
~'5'-54 

136li 
371^^ 
761^^ 


oc+260 

16211^° 

97r.+181 

9cq+126 
^00_i4i 

099+229 


288+ ^'^ 
34214^7 

410lifo 
2951 f6 
452I4II 


Qsn_4i4 

OOU 1 ± J. 

414-635 


404 
480 


29.4 
34.9 


19 
19 


129lli^ 
1559t?i 


6001^^^ 
41217/1 


12801125 
140611!^ 


63+^^ 
941^^ 


o-i 0+476 
"^^"^-337 

44311^1 


555+^^2 

30314I13 




Metal-rich sample: 117 clusters with 1.42 < (C - Ti) 


< 2.35 






35-635 


239 


17.4 


117 


133lt % 


3851 ^? 


13241 % 


7C+44 
'"-46 


1601^9^9° 


3641 


35-125 
125-207 
207-297 
297-380 
380-635 


83 
172 
258 
325 
415 


6.0 
12.5 
18.8 
23.6 
30.2 


29 
24 
22 
25 
17 


1415+111 

1199+137 

125711?^ 
149911?! 


3251 11 

Q'7Q+ 95 

344l/i^4 
304+i| 

588_276 


1379+1^^ 
1368+i^^ 

-,-174+139 
-■^-■^ '^-116 
1 971 

1357lil 


lOOl^i 
90+55 
^"J-62 

lOOl'^^ 
104l« 

901?^ 


qo + 178 
^^-256 

19311°^ 

18311?^ 
174+306 

2891^^1 


319li|'7 
345lif5 
0-1 9+ 44 

"^-"-^-138 

26814^1 
4771^?^ 


380-414 
414-635 


399 
481 


29.0 
35.0 


8 
9 


1433^1 

1512_287 


coo + 114 
000_442 

cco + 196 
000-568 


1350^ 
1350*^ 


181!^ 
116l^« 


667+ ^^"^ 
""'-1000 

200+ 


0-I4+112 
•J-L^-314 
409+36O 
'±'J^_432 



"Systemic velocity held fixed in fit of sine curve. 
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DYNAMICS OF THE M87 GLOBULAR CLUSTER SYSTEM 



of radius, particularly among the blue clusters. This interesting phenomenon is discussed in more detail below. 

3.3. Kinematic Properties as a Function of Projected Radius 

In order to better visualize these trends, and to address any concerns that our particular choice of bins in Table might 
have influenced our results, Figures show the line-of-sight mean velocity and dispersion, and the best-fit rotation 
axis and amplitude, as functions of projected radius for our full GC sample, and separately for each of the blue and red 
subsets. To construct these profiles we slide a bin of fixed radial width AR — 90" ~ 6.5 kpc through the GCS, centering 
the bin at the position of each of our measured GCs in turn and computing the kinematics (including bootstrap estimates 
of confidence intervals) for all clusters lying within ±(Ai?)/2 of that point. We begin at small projected radii when the 
number of clusters per bin first exceeds 15 (so that the fits can be usefully constrained) and we stop at large R when the 
sample size per bin drops below this limit. 

Clearly, most of the points plotted in these figures are not statistically independent — only points separated by at least 
the full AR — 6.5 kpc (the size of the horizontal bar drawn in each middle panel) can be — and we do not attempt any sort 
of formal fits to these profiles. Nevertheless, this procedure is useful since it provides a more detailed (albeit smoothed) 
view of radial trends than can be had with a single binning such as that in Table n. It is also important to note that this 
method of smoothing the data is different than that employed by Kissler-Patig & Gebhardt (1998) and Cohen (2000). 
These authors opt instead to fit the kinematics in a series of annuli that always enclose the same number of data points, 
and which therefore (because neither the intrinsic spatial distribution of the GCS nor the observational sampling of it is 
uniform) have a variable radial width (for an illustration of this point, see Figure 2 of Cohen 2000). Our approach avoids 
the nontrivial difficulty of interpreting data that have been smoothed by different amounts from point to point, and it 
yields a more accurate, though necessarily less precise, depiction of radial variations in kinematics at the low-density, 
large-radius limits of our GC sample. 

The mean velocity in Figure |^, which makes no allowance or corrections for rotation, shows no significant variations 
with radius until roughly i? > 30 kpc. The sharp apparent rise in (vp) at this radius is due — as Figure ^ and the bottom 
of Figure ||both show, and as Romanowsky & Kochanek (2000) also suggested — to sparse sampling of a population that 
shows significant rotation: the few GCs in our radial-velocity sample at these large radii have been drawn mostly from 
positions around the photometric major axis of M87 (at 155° East of North), where the rotation produces line-of-sight 
velocities furthest from the true systemic velocity of the galaxy and the GCS. Fitting a sine curve to quantify the rotation 
at every radius takes this effect into account, and a plot of Vsys vs. R (cf. eq. ^) is essentially flat, and consistent with a 
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Fig. 5. — Upper Panel: Biweight estimates for the mean line-of-sight velocity of the M87 globular cluster system, as a function of distance 
from galaxy center (filled circles). Thin and thick solid lines show the 68% and 90% confidence limits on {v), estimated using a bootstrap 
technique as described in the text. Points represent the systemic velocity measured in radial bins of constant width Afi_= 90" 2; 6.5 kpc 
(indicated by the error bar in the middle panel) . The broken horizontal line indicates the mean velocity recorded in Table hi for clusters at all 
radii 35" < R < 635". Middle Panel: Same as above, but for the sample of metal-poor (blue) clusters. Lower Panel: As above, but for the 
sample of metal-rich (red) clusters. 
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R/kpc 

Fig. 6. — Upper Panel: Biweight estimates for the velocity dispersion of the M87 globular cluster system jas a function of galactocentric 
distance. Solid points denote the dispersion Up about the sample mean in each bin, taken from Figure bl Open circles represent the 
dispersion o-p ,. about the best-fitting sine curve in each case (parameters taken from Figures M and ph. Solid curves delimit the 68% and 90% 
confidence bands for ap, obtained from a numerical bootstrap. Horizontal line is at the value determined for the full sample of clusters with 
35" < R < 635" . Middle Panel: As above, but for the sample of metal-poor (blue) globular clusters. Lower Panel: As above, for the sample 
of metal-rich (red) clusters. 



constant Vsys = 1350 km . 

The velocity dispersion profiles in Figure |^ show both the biweight rms scatter, ap, about the direct mean of the sample 
in each bin (solid points) and the biweight dispersion, crp r, about the best-fitting sine curve (open circles). At relatively 
small projected radii R ^ 10-15 kpc, the red clusters may be dynamically somewhat "hotter" than the blue clusters in 
M87, and the situation may be reversed at larger radii. As was mentioned above, however, in connection with Table 
^, these apparent differences are not statistically significant. However, in agreement with previous studies (e.g., Cohen 
& Ryzhov 1997), we do see evidence for a rising velocity dispersion towards larger radii — a trend that is more clearly 
established for our blue clusters than for the red ones. 

Figures ^ and ^ — which show the rotation-axis position angle and apparent rotation amplitude as functions of projected 
radius in the M87 CCS — present the most substantial new result of this empirical part of our analysis and constitute a 
point of clear departure from earlier studies. These profiles were obtained by fitting equation to the fp(0) data in our 
sliding radial bin, with Vsys = 1350 km s^^ held fixed but with Qq and Q.R both allowed to vary. The most striking point 
relates to the blue clusters, which appear (see Figure |^) to rotate about the projected major axis of M87 (Go = —25° 
E of N) in the inner regions of the system, R <^ 16-18 kpc, but then switch to the the more conventional minor-axis 
rotation (8o = 65° E of N) that they show for i? ^ 20 kpc. Given the smoothing kernel that we have applied to the data 
(again, AR = 90" ~ 6.5 kpc, as illustrated in the middle panel of Figure ^), this 90-degree flip in 60 seems to be almost 
instantaneous. It is probably no accident that it appears to coincide with the onset, at i? ~ 19 kpc, of the cD envelope 
in the starlight (de Vaucouleurs & Nieto 1978; Carter & Dixon 1978) and the CCS of M87 (McLaughlin, Harris, & Hanes 
1993). Figures ^ and || also reveal that the rotation axis of the red clusters seems to show a gradual drift towards the 
photometric major axis of M87 for i? ^ 18 kpc before it "corrects" to a constant Go — 65° at larger radii. It is not 
clear, however, that this apparent difference in the small-i? behavior of the red and blue GCs is statistically significant 
— even at the 68% confidence level. We return to this issue in §5 below, where the two-dimensional rotation fields of the 
metal-poor and metal-rich GC samples are discussed within the context of the larger Virgo cluster environment. 

Figure ^ shows the rotation amplitude corresponding to each fitted value of Go in Figure |^. As was suggested above, 
it is difficult to discriminate between a flat rotation curve and one which rises with increasing radius. For instance, there 
is some evidence for an outward rise in QR for i? ^ 25 kpc, for the red and blue GCs alike, but this effect sets in roughly 
where the azimuthal coverage of our radial- velocity sample begins to decline (see Figures ^ and ^ , and it is signiflcant at 
just the 68% confidence level. Also, note that the dashed curve drawn in all three panels of Figure || shows the circular 
velocity of the M87/Virgo potential in which the GCs are orbiting [vc{r) — {GM{r)/r}^/'^; see eq. ( |ll| ) below for our 
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R/kpc 

Fig. 7. — Upper Panel: Projected azimuth of the rotation axis of the full M87 globular cluster system, as a function of projected galactocentric 
radius (points). Thin and thick solid lines define GB% and 90% confidence bands. In this figure, and the following one, the systemic velocity 
was held fixed at Vsys = 1350 km s-l (see Table hi). The upper and lower horizontal lines indicate position angles of the projected major 
(©0 = —25°) and minor (©o = 65°) axes of M87. Middle Panel: For the metal-poor (blue) globular clusters. Lower Panel: For the metal-rich 
(red) globular clusters. 

adopted mass model]. The fact that flR appears to approach this hmit at our largest sampled radii suggests that the 
sparser spatial coverage there may be biasing our results to some extent; it may also provide some evidence that the M87 
GCS is being viewed nearly edge-on, as any inclination correction is likely to be small if the apparent flR is already so 
near Vc{r). 

One other point worthy of note in these fii? profiles is the apparent dip to nearly at radii R ~ 14-18 kpc. This 
may simply be an artifact of smoothing the data, rather than an indication of any real feature in the intrinsic velocity 
field of the GCS. Since 8o appears to change very suddenly by roughly 90° at some radius in this range, it is conceivable 
that the 6.5-kpc wide annular bins centered on points in this part of the GCS are simply combining two distinct samples 
of clusters, rotating about essentially orthogonal axes, whose sinusoidal Wp(6) dependences interfere destructively and 
average out to a lower amplitude of apparent rotation. 

Using the spectroscopic metallicities of Cohen, Blakeslee, & Ryzhov (1998), Kissler-Patig & Gebhardt (1998) attempted 
to distinguish between metal-poor and metal-rich clusters in their study of the rotation properties of the M87 GCS. They 
concluded that the blue clusters showed essentially no rotation, UR ~ 0, for small radii R <, 15 kpc. However, they fixed 
8o at about 30° E of N in their analysis, and solved only for the amplitude when fitting sines to their Vp{'d) data. The 
difference between their result and ours for the blue-GCS rotation amplitude stems from the fact that we have found Oq 
at small R to be far away (i.e., approaching 90°) from the position angle adopted by Kissler-Patig & Gebhardt (1998). 

The acquisition of still more radial velocities for M87 GCs will be important in confirming our findings and in clarifying 
and quantifying any dichotomies that might exist in the blue vs. red kinematics. Meanwhile, the surprisingly complicated 
rotation field that we have uncovered demands the development of sophisticated dynamical models of the GCS. In 
particular, although we proceed in the next Section with a standard, one-dimensional Jeans-equation analysis of velocity 
moments in the M87 GCS, we recognize that we cannot hope thereby to account rigorously for the dynamical effects of 
rotation. In fact, in this analysis we shall neglect rotation altogether, treating it as a perturbation to be addressed in 
future studies. 

Figure ^ roughly quantifies the effect this choice may have on our results by comparing the rotation amplitude of the 
M87 GCS to its line-of-sight velocity dispersion, both about the sample mean (filled circles) and about the best-fit sine 
(open circles), at every projected radius. The ratio of the two is essentially constant with projected radius, holding at a 
level of about 0.45 for the GCS as a whole and for the metal-poor and metal-rich subsamples individually. Thus, rotation 
enters as a correction at the [D.R/apf' ~ 20% level in the Jeans-equation analysis that we shall pursue below and, thus, 
at this exploratory stage, is a relatively small source of uncertainty. 
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4. DYNAMICAL MODELS 

In the absence of rotation and the approximation of spherical symmetry, the fundamental equation of our dynamical 
analysis is the spherical Jeans equation (e.g., Binney & Tremaine 1987): 

— TlclCT^ H n-clCTr = ~"cl 5 , (7) 

dr r 

where nci(r) is the volume density profile of the GCS, ar{r) is its velocity dispersion in the radial direction, and f3ci{r) = 
I — ag/a"^ is a measure of its velocity anisotropy.^ 

The tendency in previous studies along these lines has been to use the projected velocity dispersions of the GCs to infer 
the intrinsic ar{r) and thus to constrain the total, "background" gravitating mass profile Mtot{r). To do so, however, 
requires some a priori assumptions on the anisotropy profile, Pciif), and typically this has meant adopting the simplest 
case of isotropic orbits with /3ci = (see, e.g., Cohen & Ryzhov 1997). By contrast, Romanowsky & Kochanek (2000) 
apply a technique in which an ad hoc functional form of Mtotif) is assumed a priori, and the observed (7p{R) is used in 
conjunction with sophisticated orbit modelling to constrain both the normalization of the M87/Virgo mass profile and 
the behavior of Pciif) in the GCS. This latter procedure is closer in spirit to ours here, in which we specify both nci(r) 
and Mtot{r) in full, and then use the observed (Jp{R) profiles of §3 to infer something about the anisotropy of the GC 
orbits in M87, both globally — for our entire sample — and separately — for the metal-poor and metal-rich populations. 

One way to do this is as an 'inverse problem' (e.g. Merritt & Oh 1997). In this approach, the observed function ap{R) 
is approximated in some model-independent way (e.g., via sphnes). For an assumed Mtot(r), the two unknown functions 
o'ri^r) and f3{r) then follow uniquely from the Jeans equation and the deprojection integral. But given the limited number 
of velocities in our sample — to say nothing of our neglect of rotation and our assumption of spherical symmetry — we opt 
instead to construct a suite of model ap{R) profiles by specifying /ScK?") a priori, and then comparing these models with 
the data to gain a broad, qualitative view of the orbital properties of the M87 GCS. Formally, we solve equation (0) for 
ar{r): 



• GMtot ( [Wcx , , , 
rici 5 — exp I / dx dx 



nc\{r) \ J r 

■^^Considerations of symmetry imply that the two tangential velocity dispersions are equal, i.e., ag = cr^ 



(8) 
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Fig. 8. — Upper Panel: Amplitude of projected rotation as a function of galactocentric radius in the M87 globular cluster system (filled 
circles), shown with 68% and 90% confidence bands (thin and thick solid lines). The position angle, 0o, of the rotation axis was |aJlowed to 
vary with radius in the sinusoidal fits (see Figure M), but the systemic velocity was held fixed at Vsys = 1350 km (see Table hi). Upper 
horizontal (dash-dot) line is placed at the amplitude fiij = 169 km s~^ found by fitting to the entire sample of globular clusters, over all radii 
(see Table hi). Dotted curve shows the circular velocity as a function of radius in the M87/Virgo gravitational potential, Vc{r) = [GM{r) /r]^^'^ . 
Middle Panel: As above, but for the sample of metal-poor (blue) clusters. Lower Panel: As above, but for the sample of metal-rich (red) 
clusters. 



14 DYNAMICS OF THE M87 GLOBULAR CLUSTER SYSTEM 




10 20 30 40 

R/kpc 



Fig. 9. — Upper Panel: Ratio of the projected rotation amplitude in Figure |8| to the hne-of-sight velocity dispersion in Figure Solid 
points compare QR to the biweight dispersion dp about the mean velocity (vp) at each radius, and the solid lines trace the 68% and 90% 
confidence limits on the ratio. Open points come from comparing QR to the biweight dispersion Up about the best-fitting sine curve at each 
radius. For clarity, the confidence bands for this second ratio are not shown. Thti.upper horizontal line marks the average of ((f!/?)/(Tp) = 
0.45 suggested by a fit to all cluster in our sample, regardless of radius (see Table hi). Middle Panel: As above, but for the metal-poor (blue) 
globular clusters. Lower Panel: As above, for the metal-rich (red) globular cluster system. 



With ar(r) in hand, the projected velocity dispersion profile, ap{R), then follows from the standard integral 

where the projected GCS density Nci{R) is given by equation (|^). We also compute the cumulative, aperture velocity 
dispersion profile, defined to be the rms line-of-sight velocity for all clusters inside a given radius: 

-1 

(10) 

This statistic has the obvious advantage of reduced noise in the measured dispersion at large radii, although its interpre- 
tation can be somewhat more complicated since it reflects a coarse average (and progressively more so as R increases) of 
a spatially varying quantity. 

To proceed in this way, we require accurate descriptions of: (1) the total gravitating mass profile in M87 (and, as the 
projection integrals above make clear, in the Virgo Cluster beyond our last observed galactocentric radius); and (2) the 
spatial distribution of the full M87 GCS and of the blue and red clusters, nci{r), individually. 

First, we take Mtot{r) directly from McLaughlin (1999a), who has developed a simple mass model for M87 and Virgo 
that satisfies all observational constraints imposed by the optical surface photometry of the galaxy (as a tracer of the 
stellar mass density; de Vaucouleurs & Nieto 1978) by the X-ray surface brightness on i? ^ 200-kpc scales in the core of 
the cluster (used under the assumption that the hot intracluster gas around M87 is in hydrostatic equilibrium with the 
total mass distribution; Nulsen & Bohringer 1995), and by the spatial distribution and line-of-sight velocity dispersion 
of the early-type (mostly dE) Virgo galaxies out to scales of ~ 3 Mpc (Binggeli et al. 1985; Girardi et al. 1996). We 
emphasize that because this model for Mtot(?') was constructed explicitly without any reference to GCS data, it is 
completely separate from, and independent of, our radial-velocity observations. Thus, the total (baryonic and dark) mass 
interior to any three-dimensional radius r in Virgo is taken to be, with normalizations appropriate for a distance of 15 
Mpc to the cluster, 

Mtot(r) = M,tars(r-) +Mdark(0 

M,tar.(r) = 8.10X IQiiMo [(r/5.1kpc)/(l + r/5.1kpc)]'-^^ (11) 
MdarkW = 7.06 X 10^^ Mq [ln(l + r/560kpc) - (r/560kpc)/(l + r/560kpc)] . 



NciiR)al{R)RdR 



Nci{R)RdR 



COTE ET AL. 



15 



Note that this model attributes all of the dark matter around M87 to the Virgo Cluster as a whole; all of the available 
data are consistent with this minimalist hypothesis, in which M87 as a galaxy has no dark matter halo of its own (see 
McLaughlin 1999a). The functional form of Mdark(?') is that expected for the "universal" dark- matter halo of Navarro, 
Frenk, & White (1997), in which the volume density scales as pdark ~ f~^(''s + with a fitted value of — 560 kpc. 
The adopted models for Mstars('"), Mdavkir) and Mtot(?') are illustrated in Figure |l^. 

Next, Figure ^ motivates our favored model for the density profile of the full M87 CCS. This plot, which is taken 
nearly in full from McLaughlin (1999b) and which is described in detail there, shows the volume mass densities of the 
stars in M87 (obtained by applying a spatially constant mass-to- light ratio to the optical surface photometry of M87), 
the GCs around the galaxy (obtained by combining number counts from a number of different studies and applying 
a geometrical de-projection algorithm), and the hot, X-ray emitting gas in the core of Virgo (obtained in a model- 
independent way from ROSAT observations, by Nulsen & Bohringer 1995). The dash-dot line running through the gas 
datapoints is directly proportional to the density profile of the dark-matter halo just discussed — i.e., it shows that 
Pgas oc pdark oc (r/560 kpc) ~ ^ (1 -|- r/560kpc)~^ in Virgo (see McLaughhn 1999a). 

The GC density data in Figure ^have been scaled up by a factor of 1/0.0026 ^ 385 to illustrate that, beyond a radius 
of r « 7 kpc (the effective radius of the M87 starlight in projection) a strict proportionality ties the density of the GCS 
to that of the other baryons: Pci{r) = 0.0026 [pstars('') + Pgas{r)]. The constant of this proportionality, 0.0026, appears to 
hold similarly in nearly 100 other early-type galaxies, and it has been interpreted by McLaughlin (1999b) as a universal 
efficiency of GC formation. Thus, for radii r ^ 100 kpc, where there are no direct observations of the GCS density profile, 
we assume that Pdir) — 0.0026pgas(f ) oc pdark('')- Inside this, we adopt an ad hoc functional fit to the GCS density data 
from McLaughlin (1999b). In total, we posit that 



ntotir) 



oc 



(r/9.1kpc) 
(r/560kpc] 



r/9.1kpc) ■ 
-r/560kpc)- 



r < 95 kpc 
r > 95 kpc . 



(12) 



The normalization of this model, which is drawn as the bold, dashed line in Figure IT], is unimportant for our purposes 
(see eq. Q). The observed proportionality between the GC and gas density profiles has some important implications for 
the debate over whether some or all of the GCs surrounding M87 are intergalactic in nature (White 1987; West et al. 
1995; Harris et al. 1998). We shall return to this issue in §5.7. 

Figure ^ shows the projection of this density model as a bold solid line. (The mass surface density S is related to p in 
the same way as the projected number density N is related to n.) The small triangles are the observed surface densities, 
again taken from McLaughlin (1999b), that were de-projected to provide the volume densities in Figure pl|. We have now 
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Fig. 10. — Mass model for the central regions of the Virgo cluster, adopted from McLaughlin (1999a). The solid curve shows the total mass 
profile: i.e., the combined Virgo cluster dark matter distribution (as traced by the x-ray emitting gas) and that of the stellar component of 
M87, indicated by the dashed and dotted curves, respectively. See §4 and equation (11) for details. The vertical arrow shows the projected 
radius where the cD enevelope of M87 begins. 
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Fig. 11. — Comparison of the (mass) volume density profiles of the globular cluster system (filled circles and bold dashed line), halo stars 
(thin dashed line) and X-ray gas (open squares and dot-dash line) around M87 (after McLaughlin 1999b). Note that Pcl has been scaled 
upward by the factor 1/0.0026 (see text). The solid line is the sum of the stellar and gas densities, which have been calculated by de-projecting 
the galaxy surface photometry with the algorithm of McLaughlin (1999b). The vertical arrow marks the projected effective radius of the 
galaxy light, R^g = 96" ~ 7 kpc (de Vaucouleurs & Nieto 1978). Beyond this, it is found that Pcl ^ (Pstars + Pgas)- 

also used our new CTi photometry to measure the surface density profiles, not only for the complete GC sample, but 
also for the metal-poor and metal-rich samples separately. To do so, our initial catalog of point sources was first trimmed 
to exclude all objects having Ti < 21, C— Ti < 0.8, and C— Ti > 2.35, leaving us with a sample of 2130 candidate GCs. 
The adopted limiting magnitude was chosen to ensure that our photometric catalog is complete over the full range in 
galactocentric radius (see Geisler, Lee & Kim 2001). A background surface density of 0.8 clusters arcmin"^ was then 
subtracted to give the surface density profile for the entire cluster population. This background level was chosen based on 
a comparison of our surface density profile with the calibrated surface density profile of McLaughlin (1999b), and agrees 
to ~ 5% with that predicted by the Galactic star count model of Bahcall & Soneira (1981). The final, blue-plus-red 
profile is shown as the open circles in Figure |l2|, where it has been scaled vertically to match the counts from McLaughlin 
(1999b) and shows good agreement with that independent measurement of Sci(i?). 

For the separate metal-rich and metal-poor samples, we have taken the background surface densities to be exactly 
half that of the total background in our CTi sample. Though somewhat arbitrary, these choices of backgrounds are 
consistent with the color distribution of foreground stars predicted by Ratnatunga & Bahcall (1985) for three Galactic 
GCs having Galactic latitudes and longitudes similar to M87. This approach is obviously inferior to measuring surface 
densities directly from background fields, but the precise choice of backgrounds has relatively little effect on the final 
surface density profiles, except in the outermost bins. The resulting distributions for the metal-rich and metal-poor 
cluster populations are indicated by the open squares and circled points in Figure |l^. This shows quite clearly that 
the metal-poor clusters in M87 follow a much more shallow radial profile than their metal-rich counterparts, consistent 
with earlier findings based on HST imaging in the central regions of the galaxy (e.g., Neilsen, Tsvetanov & Ford 1999). 
There is no reasonable choice of background that can alter this basic conclusion. Even if we were to make the highly 
implausible assumption that the background consists entirely of sources having colors of C— Ti < 1.42, this qualitative 
contrast between the blue and red GCSs would remain. 

A more detailed discussion of the spatial distribution of the M87 GCS will be presented in Geisler, Lee & Kim (2001); 
our purposes only require estimations of the overall shapes of rtbiuo('") and Tired (?")• To obtain these, we project models 
similar in form to that in equation (03) and fit to the blue- and red-GCS data points in Figure O. We obtain 
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cx (r/560kpc)-i(l + r/560kpc)~2 ^ r > 95 kpc , 



(13) 



COTE ET AL. 



17 




I I I I 1 1 I I I I I I I I I 

1 10 100 

R/kpc 

Fig. 12. — Surface density profiles, for globular clusters associated with M87. Triangles show the surface density profile of McLaughlin 

(1999b) for the full globular cluster system, while open circles show the profile found here using the Washington photometry described in 
Paper I. Circled points and open squares indicate the respective surface density profiles of metal-poor (blue) and metal-rich (red) globular 
clusters. The bold curves correspond to projections of the models in equations (12)— (14) in the text; see the discussion there for details. 



and 

nred{r) oc (r/3.3kpc)-i(l + r/3.3kpc)-2 , r < 125 kpc , - 

(X (r/560kpc)-i(l-f r/560kpc)-2 , r > 125 kpc . ^ ' 

The projections of these best-fit models are drawn as the bold, broken lines in Figure ll^. Their absolute normalization 
is again unimportant for us, although their relative normalization does matter. This hasbeen set so that rtbiuc/'T-rcd 4, 
independent of radius in the limit r — s- oo. The projected sum of these red and blue models is shown in Figure ^ as the 
thin solid line, which is reasonably compatible with the completely independent projection of n^ati^) from equation (E2[). 

With specifications for the Virgo mass profile Mtot{i') and for the density profiles rici(r) of the total, blue, and red GCSs 
thus in hand, we solve equations (P)-(]l0|) for five different spatially constant values of the GCS orbital anisotropy Pdf^ 



-0.99 






0.01, radial bias) 


+0.4 




lol^ 


0.6, radial bias) 





(erg 


~ cr^ 


— (Tr, isotropic) 


-0.4 




lol^ 


1.4, tangential bias) 


-99 




lol = 


100, tangential bias) 



/3ci = < {oq = 0-0 = (Tr, isotropic) (15) 



Figures 13 -IS] show the results of this exercise applied to our three GCS samples. In each figure, the top panel shows 
the observed <7p(K) profile computed by smoothing the data with a radial kernel of fixed width 120" ~ 8.2 kpc (rather 
than the 90" used for Figure ^). Also shown are the 90% confidence bands, obtained by bootstrapping as described in §3. 
The bottom panels show the aperture velocity dispersion and bootstrap estimates of the 90% confidence bands for the 
cluster color ranges indicated. Overlaid on the data in every case are the model curves predicted by substituting equations 
([ll|)-(p^) in equations ([8|)-([To|). We emphasize that these theoretical velocity dispersion profiles have not been fit to the 
data: with Mtot(?'), ?^ci('^), and /3ci all specified a priori, there are no additional parameters to be adjusted. Note that 
the rising velocity dispersion profile seen in Figure ^ is once again evident here — not only in the measured ap (R) and 
o'ap(^) profiles, but also in the models. This predicted rise in a^R) and crap(^) comes about because the background 
mass density profile given by equation |ll], and shown in Figure |lO|, is shallower at these radii than that of an isothermal 
sphere. 

It is clear from Figure ^ that the M87 GCS as a whole has an almost perfectly isotropic velocity ellipsoid. Models with 
even modest amounts of anisotropy. Pel = ±0.4, produce noticeably inferior matches to observed velocity dispersions, while 

^^It is worth emphasizing that these values of /S^i are best interpreted as density-weighted averages due to our assumption of spatially constant 
anisotropies. 
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Fig. 13. — Upper Panel: Velocity dispersion profile for the complete sample of globular clusters. Points represent biweight estimates for 
the line-of-sight velocity dispersion in radial bins of width AR = 120" ~ 8.7 kpc; thin solid curves show bootstrap estimates of the 90% 
confidence limits on these measurements. The thick solid line shows the predicted velocity dispersion profile for a system of test particles that 
are embedded in the potential defined by equation (11), follow a density profile given by equation (12) (the bold, dashed line in Figure |l]| ), 
and have an isotropic velocity ellipsoid (/3(.i = 0). From top to bottom, the four dashed curves show the predicted velocity dispersion pronles 
for = +0.99 (a strong radial bias), /3;.i = +0.4 (a moderate radial bias), /3^i = —0.4 (a moderate tangential bias) and /3^i = —99 (a strong 
tangential bias). Lower Panel: Aperture velocity dispersion profile for the complete sample of globular clusters, with bootstrap estimates of 
the 90% confidence limits illustrated by thin solid lines. The heavy solid curve shows the predicted aperture velocity dispersion profiles for 
the case of isotropic orbits, while the four dashed curves show model predictions with the same velocity anisotropies as in the above panel. 



strongly tangentially biased orbits (/3ci — —99) can be ruled out with > 90% confidence and very radial orbits (/3ci — 0.99) 
are altogether out of the question. However, if the cluster sample is divided on the basis of metallicity, some interesting 
differences emerge. Figure O suggests a modest tangential bias in the metal-poor clusters, with (3ci ~ —0.4 at small radii 
(corresponding to ae ~ 1.2ov), perhaps tending towards isotropy at larger radii '^20-30 kpc. The opposite result is seen 
for the metal-rich GCs, which have a slight radial bias of roughly the same magnitude: Figure |l^ suggests f3ci ~ -1-0.4, 
and thus ag ~ O.Sar- This is just what is needed, of course, to balance the anisotropy in the blue GCS and produce the 
isotropy indicated for the combined sample in Figure |l^. Qualitatively speaking, these contrasting anisotropies are also 
what might have been expected a priori, given the fact that the blue GCS has a much shallower density profile than the 
red subsystem, even though the two populations are embedded in the same gravitational potential and show no significant 
differences in their projected kinematics (see Figure 

Romanowsky & Kochanek (2000) have derived a (3c\{r) profile for the total GCS — i.e., they did not separate their 
velocity sample on the basis of metallicity. At small radii, the implied (3ci is markedly negative, thereby implying a 
substantial tangential bias for the GCS. However, this claim ultimately rests on their assumption that the background 
potential of M87/ Virgo is distributed as a singular isothermal sphere, Mtot{r) oc r, whereas the true Mtot{r) in the core 
of Virgo is well constrained (principally by X-ray data; see Nulsen & Bohringer 1995 and McLaughlin 1999a; see also 
Figure ^0|) to grow substantially more rapidly than this with increasing radius. Said another way, the density distribution 
of the gravitating mass in the core of Virgo is significantly shallower than the p cx distribution of an isothermal sphere, 
and therefore the density distribution of the central GCS does not require so large a tangential bias as Romanowsky & 
Kochanek (2000) found necessary to satisfy the Jeans equation. 

Finally, we consider the effect on our results if our adopted extrapolation of the GCS density profiles in Figure n3 and 
equations (|l^-(0) is somehow in error. It seems unlikely that the GC distribution could become shallower than me hot 
intracluster gas (and the dark matter) in Virgo, and thus Figures |l6| and ^ explore the implications of an nci(r) that 
drops off very steeply towards large radii. Rather than requiring n^r) oc Pgasir) at the unobserved r ^ 100 kpc, then, 
we take at face value simple fits to our CTi density profiles for 5 ^ i? ^ 30 kpc and assume 



nuucir) oc (r/21.3kpc)"i(l + r/21.3kpc)"2 , aU radn 
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Fig. 14. — Same as Figure hA except for the sample of metal-poor (blue) globular clusters. Note that the aperture velocity dispersion at 
small radii in the lower panel rails below that expected in the case of isotropic orbits, at the 90% confidence level. The model which best 
matches the measured profile has /B^i —0.4, corresponding to a tangential bias of ag fsi 1.2o'r- 
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Fig. 15. — Same as Figure [l3| except for the sample of metal-rich (red) globular clusters. In this case, the aperture velocity dispersion 
profile at large radii in the lower panel lies above that expected for an isotropic velocity ellipsoid. A model which reasonably matches the 
measured profile has f3^i « +0.4, corresponding to a radial bias of ag 0.8(Tr. 
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Fig. _L6. — Alternate extrapolation of the red and blue globular cluster system density profiles to larger, unobserved radii (compare with 
Figure hd). Thin, solid line is the sum of the red and blue fits (shown as thick lines) and is in good agreement with the independently observed 
density profile of the full globular cluster system. 



and 

rircdir) oc (r/3.3kpc)"^(l + r/3.3kpc)"2 , all radii. (17) 

The bold, broken lines in Figure |l^ compares the projection of these functions to the same data considered in Figure ^ 
The thin, solid line, which shows the sum of the red and blue projected densities, is again in good agreement with the 
independently observed density distribution of the full GCS. 

We use these new density profiles with the same M^ot (f) from equation ( pT| ) and the same five constant (3^1 values from 
equation jisl ) to compute new models for the (Tp{R) and (Tap(< R) of the metal-poor and the metal-rich components of 
the GCS. Figure |l^ compares these to the data, which again have been smoothed with a 120"-wide sliding radial bin and 
are shown with bootstrap estimates of their 90% confidence bands. Because the large-r density gradient of both cluster 
subsystems has been made steeper by hypothesis without changing any other aspect of the dynamical models, the data 
would appear now to suggest a move towards substantial radial bias in orbits of both blue and red clusters at the largest 
radii covered here. We stress, however, that we regard this particular scenario as a very extreme possibility, as the GCS 
density profiles postulated earlier are certainly closer to the truth. For instance, if nci(r) indeed behaves as in Figure |l^ 
rather than as in Figure ^ then we must conclude that the GC-to-baryon mass ratio of 0.0026 is only coincidentally 
realized, and then over a limited range of radius, in M87 — even though this "universal" ratio matches precisely what 
is found by McLaughlin (1999b) for nearly 100 other early-type galaxies. In any case. Figure |l^ also shows that the 
GC volume density profiles cannot fall off more rapidly than roughly nd oc r~^ , since this gradient already requires the 
clusters at large r to have nearly purely radial velocity ellipsoids in order to reproduce the observed line-of-sight velocity 
dispersion. 

5. DISCUSSION 

5.1. Results for other Ellipticals from Globular Clusters 

Before turning our attention to the various models which have been proposed for the formation of GCSs and their host 
galaxies, we pause first to compare our findings for the GCS of M87 with those of other elliptical galaxies. Perhaps the 
most natural comparison is that of M49, the most luminous member of the Virgo cluster, whose GCS has been studied 
recently by Sharpies et al. (1998) and Zepf et al. (2000). From a sample of 144 GCs with measured velocities, Zepf 
et al. (2000) report that the metal-poor and metal-rich GCSs have different dynamical properties: i.e., frot/f ~ 0.3 and 
0.1, respectively (with Wrot/c < 0.56 and 0.34 at 99% confidence). These estimates are to be compared with the value of 
t^rot/c ~ 0.45 which we measure for both components in M87. Hui et al. (1995) have analysed the radial velocities of 
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433 planetary nebulae (PNe) and 62 globular clusters in NGC 5128, and conclude that the metal-rich GCs in that galaxy 
show roughly the same dynamical properties as do the PNe, Wrot/o' ^ 0.3, whereas the metal-poor GCs show no evidence 
of rotation. Meanwhile, Kissler-Patig et al. (1999) find some marginal evidence for rotation among the outermost GCs 
in NGC 1399, the central elliptical in Fornax, but due to the limited size of their radial velocity sample, were unable 
to ascertain if the chemically distinct GC populations differ in their dynamical properties. Thus, at this early stage, no 
obvious pattern in the dynamical properties of the GCSs of elliptical galaxies has yet emerged. 

5.2. Results for other Ellipticals from Planetary Nebulae 

Radial velocity surveys of PNe represent an alternative, and complementary, method of studying the halo dynamics of 
giant elliptical galaxies. While there have been no dynamical studies of the PNe population of M87 itself, some comparisons 
are possible for the three luminous early-type galaxies having large numbers of measured PNe radial velocities: NGC 5128 
(Hui et al. 1995), NGC 4406 (Arnaboldi et al. 1996) and NGC 1316 (Arnaboldi et al. 1998). The PNe in all three galaxies 
show evidence for significant rotation, Wrot/c ^ 0.5-1, in their outer regions. In this respect, these galaxies resemble M87 
which has {{^lR)/a) ~ 0.5 from its GCS, with some evidence for slightly lower values in the inner regions of the galaxy 
and an outward increase in {VLR)/a. Moreover, integrated-light spectroscopy in the inner regions of these galaxies reveals 
little or no rotation in each case, similar to the situation for luminous elliptical galaxies, and for M87 in particular (e.g., 
Vrot,*/<^* ^ 0.1 according to Davies & Birkinshaw 1988). On the other hand, none of these galaxies has the generally 
round and smooth photometric appearance of M87: NGC 5128 and NGC 1316 are clear examples of recent gas-rich 
mergers or accretions (e.g., Schweizer 1980), and NGC 4406 is probably best classified as an SO galaxy (Arnaboldi et 
al. 1996). While firm conclusions must await the measurement of radial velocities for PNe in the vicinity of M87, our 
results for the M87 GCS appear consistent with the claim of Arnaboldi et al. (1998) that giant elliptical galaxies, by 
virtue of their relatively rapid rotation at large radii, may contain roughly as much angular momentum per unit mass as 
do spiral galaxies assuming that the GC and PNe kinematics are representative of the kinematics of the underlying galaxy 
mass. 

5.3. The GCS Velocity Ellipsoid at Small Galactocentric Radii 

It has been known for some time that the surface density profile of the GCs surrounding M87 exhibits a "core" of 
radius Rc — 4-5 kpc (Lauer & Kormendy 1985; McLaughlin 1995). The possibility that this core developed as a result 
of GC destruction through tidal shocks from a central bulge or black-hole, or from an underlying triaxial potential, has 
been explored by many investigators (e.g., Ostriker, Binney & Saha 1989; Pesce, Capuzzo-Dolcetta & Vietri 1992). These 
authors also noted that such processes are expected to lead to a preponderance of tube orbits among the surviving GCs 
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Fig. 17. — Model vs. observed dispersion pro files shown separately for the blue and red subsets of the M87 globular cluster system, assuming 
the alternate density profiles shown in Figure llQ and described by equations (16) and (17). Since the mass profile of M87/Virgo is fixed in 
this analysis, the data in this case prefer modeTs with stronger radial bias in the globular cluster system orbits at large galactocentric radii. 
Curves shown have constant f} = 0.99, 0.4, 0, —0.4 and —99, as in Figures |l3|-|l5[) 
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relative to the initial sample. Moreover, in the specific case of the Galactic GCS, Murali & Weinberg (1997) have shown 
that GCs on eccentric orbits will undergo the most rapid relaxation and evaporation due to enhanced tidal heating; the 
preferential destruction of such clusters will therefore lead to an orbital distribution which becomes more tangentially 
biased with time. Thus, depending of course on the initial shape of the GCS velocity ellipsoid, the tangential bias which 



we find for the GCs at small galactocentric radii — especially within an effective radius (see Figures 13-151) — may provide 
evidence that these dynamical processes have played a role in shaping the presently observed GCS. 

5.4. Cooling Flow Models 

Fabian, Nulsen & Canizares (1984) were the first to propose that the populous GCSs associated with some cD galaxies 
might be the signature of GC formation through cooling flows from their X-ray halos. A number of difficulties with this 
scenario have been raised over the years: e.g., most "high-SAr" and/or cooling- flow galaxies show no evidence for young 
clusters, the GCs are generally more metal-deficient than the intracluster gas, and it is unclear how the cooling gas might 
condense into clumps of mass ^ 10^ Mq. For M87, which is both an incipient cD galaxy and the prototypical "high-SAr" 
system, we may now add another difficulty to this list: as described in §4, the orbits of the surrounding GCs are closely 
isotropic, not the radial orbits expected if the GCs have condensed out of the inf ailing gas. As Figures |l3|-|l5| demonstrate, 
purely radial orbits are ruled out at very high significance by our observations. Thus, it appears highly unlikely that a 
large fraction of M87's GCs could have formed by mass drop-out from a cooling flow. 

5.5. Monolithic Collapse Models 

Can the relatively rapid rotation of the M87 GCS be accommodated within the context of formation models which 
argue for the monolithic collapse of a single protogalactic gas cloud? Like those of most other giant elliptical galaxies, 
the M87 GCS is disjoint in terms of its chemical properties (i.e., the GC metallicity distribution function is decidedly 
bimodal; Whitmore et al. 1995; Geisler et al. 2001; see also Figure |l|), and this observation alone provides a compelling 
argument against the monolithic collapse picture. Figures may provide additional evidence against this model: both 
the metal-poor and metal-rich GC subsamples exhibit significant rotation, with QR ~ 170 km s~^ and {{flR)/ap) = 0.45. 
In monolithic collapse scenarios, the angular momentum of galaxies can only arise from tidal torques from companions 
(Peebles 1969). A common measure of a galaxy's angular momentum content is the dimensionless spin parameter, 
A = J\E\^/^G-^M-^/^, where J, E and M are the angular momentum, binding energy and mass of the collapsed galaxy. 
Fall (1979) showed that, for an elliptical galaxy in gravitational equilibrium, A ~ 0.3((fJi?)/crp). Thus, the inferred spin 
parameter of the M87 GCS is A ~ 0.15 — lower than, but still consistent with, the value of A ~ 0.18 found by Kissler-Patig 
& Gebhardt (1998) from a smaller sample of GCs. This spin parameter is roughly twice that expected from analytical and 
numerical simulations of the collapse of a single, isolated protogalactic cloud (Peebles 1969; Efstathiou & Jones 1979). It 
is, however, worth bearing in mind that M87 is by no means an isolated system. Perhaps a more severe obstacle for such 
models is the abrupt shift in the orientation of the rotation axis for the metal-poor GC subsample at i? ~ 16-18 kpc. If 
real, this feature would be difficult to explain in models where M87 and its GCS formed via the collapse and spin-up of a 
single protogalactic cloud. 

5.6. Merger and Accretion Models 

Two variations on the general theme of mergers and accretions have been proposed for the formation of globular cluster 
systems: (1) major mergers involving a pair of gas-rich disk galaxies (each containing only metal-poor GCs) which results 
in the formation of new, metal-rich GCs (Ashman & Zepf 1992); and (2) hierarchical growth of a pre-existing "seed" 
galaxy (containing its own metal-rich GCS) through the accretion of numerous smaller galaxies or protogalactic fragments 
which harbor predominantly metal-poor GCs (Cote, Marzke & West 1998). We discuss the viability of these two models 
in light of the dynamical evidence presented above. 

We begin by reiterating the well established result that mergers of purely stellar disks are a highly unlikely mechanism 
for the production of giant elliptical galaxies. Both N-body simulations {e.g., Barnes 1992; Hernquist 1992) and simple 
arguments regarding central phase-space densities (Carlberg 1986) have demonstrated that the end-products of such 
mergers have cores which are far more diffuse than the centers of giant elliptical galaxies. In addition, luminous elliptical 
galaxies typically have Hubble Types of E2 (Franx, lUingworth & de Zeeuw 1991), whereas the end-products of pairs of 
stellar disks are frequently more elongated, with Hubble Types of E3-E7 (Hernquist 1992)]^ Finally, the merger end- 
products often show large misalignments between their rotation and minor axes (Barnes 1992; Weil & Hernquist 1994; 
Heyl, Hernquist & Spergel 1996; Weil & Hernquist 1996). From an observational perspective, such large misalignments are 
rare: Franx, lUingworth & de Zeeuw (1991) find much smaller misalignments for the majority of giant elliptical galaxies 
based on long-slit spectroscopy of their integrated light. 

There are various ways to surmount some or all of these difficulties within the context of the merger hypothesis: (1) 
the stellar disks may initially harbor compact bulges with high phase-space densities (Hernquist 1993); (2) the disks 
may contain significant amounts of gas (Barnes & Hernquist 1996); or (3) ellipticals may be the result of multiple, 
dissipationless mergers (Weil & Hernquist 1994; 1996). From a purely philosophical perspective, the first option is 
obviously unsatisfactory given that the initial conditions for spiral-spiral mergers require a significant fraction of the 
initial mass to reside in compact, dynamically hot components. In the words of Hernquist (1993), "It seems likely that 

^•'And M87 itself has a Hubble^ Type of E2 (Van der Marel 1991), though the ellipticity does increase outward, reaching a maximum of E4 
in the cD envelope (see Figure hq) . 
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even the most strident critics of the merger hypothesis would admit the possibihty that mergers of small ellipticals will 
yield big ellipticals." 

The second alternative — that giant ellipticals form during dissipative mergers of pairs of stellar disks, each of which 
contain significant gas components — is closely related to the Ashman & Zepf (1992) model for the formation of GCSs 
in giant elliptical galaxies. In the specific case of M87, the most compelling piece of evidence for a past major merger 
is the possible flip in the position angle of the rotation axis for the metal-poor GCS; we caution, however, that this 
interpretation of the observed discontinuity is not unique (see §5.7). As a further complication, the uncertainties involved 
in correctly modeling the effects of gas cooling, feedback and star formation in dissipational spiral-spiral simulations are 
formidable. To date, the most detailed study of merging galaxies to include gas-dynamical effects is that of Barnes & 
Hernquist (1996), who approximated the interstellar medium as an incompressible fluid. They found that the inclusion 
of a gaseous component has a dramatic effect on the properties of the merger end-products, and may alleviate some of 
the problems described above for the case of purely stellar mergers: i.e., including gas leads to merger remnants which 
have somewhat higher central densities, and may also produce closer alignments between the angular momenta and minor 
axes. These conclusions, however, are highly dependent on how the gas is assumed to cool (see, e.g., Mihos & Hernquist 
1994, and especially §4.3 of Barnes & Hernquist 1996). 

Simulations of merging stellar disks have shown that the angular momentum axes of the merger remnants often show 
large misalignments with the minor axes. At large radii, metal-poor and metal-rich GCSs in M87 both rotate around axes 
which are closely aligned with the minor axis. If the GCs are accurately tracing the kinematics of the underlying galaxy 
mass, this alignment may constitute a difficulty for the major merger model; that is to say, such a small misalignment in 
the merger remnant is by no means implausible in the major merger model, but neither is it predicted. In their simulations 
of gas-free mergers of equal-mass disks, Heyl, Hernquist & Spergel (1996) find the distribution of misalignment angles 
to be nearly flat, with occasional misalignments of 75° or 80° although, as mentioned above, gas may play a role in 
ameliorating these large misalignments. 

It is clear that additional research into the influence of gas on the non-axisymmetry of merger remnants is needed 
urgently. Moreover, it is important to quantify the expected angular momentum content and rotational properties of 
any stars or GCs which might form during dissipational mergers. In the case of M87, both the metal-poor and metal- 
rich samples show significant rotation, with {{VlR)/ap) ^ 0.45. This value is intermediate to the disparate values of 
{{nR)/ap) - 1 (Cote 1999) and {{Q.R)/ap) - 0.1 (Zepf et al. 2000) observed for the metal-rich GCs belonging to the 
Milky Way and NGC 4472, respectively. Yet the metal-rich GCs in both of these galaxies have been identified as the 
end-products of gas-rich mergers (i.e.. Ashman & Zepf 1992 and Zepf & Ashman 1993 for the Milky Way; and Zepf et al. 
2000 for NGC 4472). 

In the dissipationless hierarchical growth scenario, the simulations of Weil & Hernquist (1994; 1996) indicate that 
the rotation axis of the merger remnant is expected to be more closely aligned with the photometric minor axis.Q 
Strictly speaking, these simulations — which do not include gas — apply to merger remnants of small virialized groups 
of galaxies, although the generic kinematic properties of the end-products are similar to those found in cosmological 
simulations of galaxy formation through the hierarchical agglomeration of smaller galaxies and protogalactic fragments 
(Barnes & Efstathiou 1987; Frenk et al. 1988; Quinn & Zurek 1988; Warren et al. 1992; c.f. Dubinski 1998). In general, 
the assumed initial conditions would be appropriate if centrally dominant galaxies such as M87 are the result of mergers 
of small subgroups which later coalesced into larger structures, as originally suggested by White (1982). The end-products 
of these simulations appear nearly round in most projections, and exhibit small kinematic misalignments when integrated 
over the full range in radius. Near their centers, they show relatively little rotation, but can rotate rapidly in their outer 
regions, with <^ t'rot/f < 0.8, although this conclusion is likely to depend on the assumed orbital properties of the 
progenitor galaxies (Dubinski 1998). 

In short, these results are broadly consistent with the observed behavior of the metal-poor GCS in M87, but the 
poorly-constrained rotation curve at large radii makes definite conclusions impossible. There is also some evidence that 
the inferred rotation of the M87 GCS (particularly that of the metal-poor component) may include a contribution from 
large-scale streaming motions along the Virgo principal axis (see below). And while the simulations of Weil & Hernquist 
(1996) produce end-products with low triaxialities and roughly oblate shapes (consistent with the results of dynamical 
modeling of the line-of-sight velocity profile in the inner regions of M87; e.g., van der Marel 1991; 1994; Bender et al. 
1994), it is important to bear in mind that the N-body simulations, though similar in some respects to the hierarchical 
growth model of Cote et al. (1998), are not completely analogous. Most notably, the N-body simulations do not include 
a single dominant progenitor galaxy, whereas in the agglomeration picture of Cote et al. (1998), the metal-rich clusters 
arise in the dissipative (monolithic) collapse of this component. 

A more definitive test of these scenarios will require a crucial piece of evidence which is currently lacking: accurate age 
estimates for the metal-rich and metal-poor GCs. To date, the limited observational evidence have proven inconclusive. 
Kundu et al. (1998) argued, on the basis of broadband VI photometry, that the metal-rich GCs in M87 are 3-6 Gyr 
younger than their metal-poor counterparts, as expected in the major merger model model. However, this claimed age 
difference is roughly the same size as the uncertainties involved in using broadband VI colors to derive GC ages, and is 
in contradiction with the findings of Cohen et al. (1998) who, from Keck spectroscopy of the brightest clusters in M87, 
found the metal-rich and metal-poor GC populations to be both old and coeval (i.e., T 13 Gyr). 

■^^Direct evidence that this process has played a role in the formation and evolution of M87 has been presented by Weil, Bland-Hawthorn & 
Malin (1997) who identified a diffuse plume of stellar material at large radii — the likely result of the accretion of a small satellite galaxy. 



24 



DYNAMICS OF THE M87 GLOBULAR CLUSTER SYSTEM 



The principal obstacle to understanding the merger/ accretion history of M87 remains the uncertain role played by gas. 
From a theoretical point of view, improved rnimerical simulations that include gas-dynamical effects and star formation 
are needed urgently. On the observational front, it is important to measure the star formation history of the GCS directly, 
through improved age estimates for the chemically distinct GC populations (e.g., Jordan et al. 2002). 

5.7. Connection to the Virgo Cluster on Larger Scales 

It has long been known that the brightest of Virgo's giant elliptical galaxies define a remarkably linear pattern on the 

plane of the sky (Arp 1968). It was recently shown that this extended arrangement, which passes through the cluster core 
at a position angle of 110-125° and is roughly centered on M87, is the projection of a highly coUinear, three-dimensional 
filament - the so-called "principal axis" of the Virgo cluster (West & Blakeslee 2000). Thus, the rotation axis of the M87 
globular cluster system not only coincides with the galaxy's minor axis at large galactocentric radii, but is also roughly 
aligned with the minor axis of the surrounding cluster. It is then natural to ask if some fraction of the GCS around M87 
may somehow be associated with Virgo itself (White 1987; West et al. 1995). One mechanism that might give rise to 
such a population of "intracluster" GCs is their acquisition from other galaxies by tidal stripping (White 1987; West et al. 
1995). Another is the possible formation of GCs in "failed dwarfs" throughout Virgo: low-mass, loosely bound clouds of 
gas — analogous to, but rather more isolated than, the protogalactic fragments thought to form stars and globular clusters 
in most largo galaxies (see, e.g., Harris & Pudritz 1994; Cote et al. 1998, 2000) — that formed normal complements of 
(generally metal-poor) globular clusters, but were then destroyed before much further star formation occurred, perhaps 
as a result of galactic winds driven at least in part by the GC formation itself (Harris et al. 1998; McLaughlin 2000). 
Whether or not the numbers of any such clouds were enhanced in Virgo through cold-dark-matter biasing of either the 
GCs (West 1993) or their host protogalactic fragments (see Ferguson & Binggeli 1994), the plausible result could have 
been an addition to the intracluster medium of globular clusters and hot gas that eventually settled naturally to the 
bottom of the cluster's potential well, i.e., directly around M87. 

The evidence in support of such a connection between the M87 GCS and the Virgo cluster can be summarized as 
follows. First, like all cD galaxies, M87 possesses an extended, low-surface-brightness envelope, and much evidence exists 
that such envelopes are entities quite distinct from cD galaxy cores. For instance, these galaxies obey the Dn-cr relation 
defined by normal ellipticals only after the luminosity excess above an i?^/^-law fit to the surface brightness profile of the 
interior is excluded (Tonry 1987); the total luminosity of cD envelopes correlates with several cluster properties including 
cluster richness and X-ray luminosity; the surface brightness profiles of the envelopes obey density profiles that roughly 
match that of the galaxies in their parent clusters (Schombert 1988); and the ellipticity and position angle of the galaxy 
isophotcs often show dramatic discontinuities at the radius where the envelope begins to dominate the siirfacc brightness 
profile. This latter issue is particularly germane in the case of M87. McLaughlin, Harris & Hanes (1993, 1994) showed that 



Fig. 18. — Orientation of the angulcir momenta axes of metal-poor and metal-rich globular clusters surrounding M87 (blue and red arrows, 
respectively) overlaid on an image from the Palomar Digital Sky Survey. A fixed radial bin width of AR = 90" ~ 6.5 kpc has been adopted 
throughout. The length of each arrow is proportional to the amplitude of the best-fit sine curve in each annulus. The maximum positive 
rotation velocity at each radius is found 90° counter-clockwise of the direction indicated by the arrows. The image measures 20' x 20' on a 
side, or 87.3x87.3 kpc for our adopted M87 distance. The luminosity-weighted mean photometric minor and major axes for the galaxy are 
indicated by the solid yellow lines. The dotted yellow curves indicate the best-fit ellipses to the surface brightness profile from Carter and 
Dixon (1978). The dashed ellipse indicates the approximate location and orientation of the galaxy's cD envelope. 
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Fig. 19. — Panel a: The line-of-sight, rotational velocity field for tlie M87 globular cluster system obtained using the full sample of 278 
clusters with measured radial velocities (dots). North is up and east is to the left. Contours are labeled in km s~^ with respect to the mean 
velocity. Solid contours indicate positive velocities, while dashed contours indicate negative velocities. Panel b: Same as (a) except for the 
sample of 117 metal-rich (red) globular clusters. Panel c: Same as (a) except for the sample of 161 metal-poor (blue) globular clusters. Panel 
d: Plane-of-best-fit to the rotational velocity field of metal-poor globular clusters. Contours are labelled in km s~^ with respect to the mean 
velocity. As in (a), solid contours indicate positive velocities, while dashed contours indicate negative velocities. Panel e: The line-of-sight, 
rotational velocity field for 117 Virgo galaxies located within 2° of M87 (dots). North is up and east is to the left. Contours are labelled in km 
s~^ with respect to the mean velocity. As in (a), solid contours indicate positive velocities, while dashed contours indicate negative velocities. 
Panel f: Plane-of-best-fit to the rotational velocity field of Virgo galaxies located within 2° of M87. Contours are labelled in km s~^ with 
respect to the mean velocity. As in (a), solid contours indicate positive velocities, while dashed contours indicate negative velocities. 



the GCS itself exhibits a sharp jump in surface density which coincides with the onset of the stellar cD envelope, and that 
the ellipticity and orientation of the GCS roughly matches that of the underlying envelope. Moreover, as was described 
in §4, the so-called "excess" of globular clusters at large radii in M87 is now recognized to be not a genuine surfeit, but a 
consequence of the fact that the GCS-to-baryon mass ratio there has the same value, Pci{r)/[pgas{i") + Pstars('')] = 0.0026, 
that is exhibited by nearly one hundred other gala xies (McLaughlin 1999b). Since the stellar component of M87 dies out 
relatively rapidly beyond r ^ 50 kpc (e.g., Fig. |ll|), the natural interpretation is that the distant GCs are associated with 
the X-ray gas around M87 — and that gas is itself at the virial temperature of the Virgo cluster dark matter halo (see, 
e.g., McLaughlin 1999a). 

The dynamical analysis presented in this paper adds to this line of argument. In particular, it is noteworthy that the 
position angle of the rotation axis of the metal-poor GCs shows an apparent discontinuity close to the onset of the cD 
envelope. To aid in the visualization of this discontinuity. Figure 18 shows the orientation of the angular momenta axes 
of metal-poor and metal-rich GCs surrounding M87 overlaid on an image from the Palomar Digital Sky Survey. The blue 
and red arrows refer to the metal-poor and metal- rich GCs, respectively. The length of each arrow is proportional to the 



26 DYNAMICS OF THE M87 GLOBULAR CLUSTER SYSTEM 



1000 




10 100 1000 

R/kpc 



Fig. 20. — Aperture dispersion profiles of the metal-poor (blue) and metal-rich (red) globular clusters around M87, and of the early-type 
galaxies in the Virgo Cluster at large radii (data from Girardi et al. 1996, shown with 95% confidence limits as erronJaars). Curves are the 
full extent of the models for the blue glohuJ^r cluster system with a fixed orbital anisotropy of /3 = —0.4 (cf. Figurea.h4), for the red globular 
cluster system with /3 = -1-0.4 (cf. Figure hsh, and for the two subsamples combined with a fixed /3 = (cf. Figure haj"^ 



amplitude {flR) of the best-fit sine curve in radial bins of fixed width AR = 90" ~ 6.5 kpc (see Fig. H). The dotted yellow 
curves show the best-fit ellipses to the surface brightness profile of the galaxy, taken from Carter and Dixon (1978). The 
ellipse nearest the onset of the cD envelope is shown by the bold, dashed curve. Recall from §4 that in the limit of r ^ oo, 
'^biuc/'^red ~^ 4. Thus, the metal-poor component dominates the GCS at these radii, although both components (i.e., the 
metal-rich and the metal-poor cluster systems) appear to rotate around an axis which is roughly aligned with the minor 
axis of the large-scale structure in Virgo. At the largest observed radii, the rotation curve of the GCS approaches the 
circular velocity of Virgo itself (although the uncertainties in the fitted values of flR are large at these distances) . 

To examine the apparent rotation of the M87 GCS in more detail, and to explore the possible connections with bulk 
motions in the Virgo cluster, we have generated two-dimensional maps of the GC velocity fields using the non-parametric 
techniques described in Gebhardt et al. (1995) and Merritt, Meylan & Mayor (1997). In this approach, the cluster-by- 
cluster velocity residuals about a mean, line-of-sight velocity are smoothed with a "thin-plate smoothing spline" which is 
chosen on the basis of a generalized cross validation technique (Wahba 1990). The results of this exercise are presented in 
Figure |l^. Panels (a-c) in this figure show smoothed velocity fields for the full sample of GCs, the metal-rich subsample 
and the metal-poor subsample, respectively. A striking feature of this figure (panel b) is the form of the rotational field for 
the metal-rich GCs: i.e., it appears non-cylindrical in nature, consistent with the discussion around equations (^-(||) in §3. 
Rather, the rotation field shows a "double-lobed" pattern with maxima at i? ~ 3.5-4i?e (25-30 kpc) along the approximate 
photometric major axis of the galaxy, although additional radial velocities may be needed to fully characterize the true 
two-dimensional rotation field. Panel (c) suggests that the rotation field of the metal-poor GCS may more closely resemble 
that of a solid body; or, alternatively, it may show evidence for a "shear" in the line-of-sight velocity. The magnitude and 
orientation of this putative shear is apparent in panel (d) of this Figure, which shows the plane of best-fit to the velocity 
residuals for the metal-poor GCS. Interestingly, the sharp flip in the rotation axis of the metal-poor clusters seen in the 
previous Figure is not obvious in this (smoothed) representation of the two-dimensional velocity field, although there does 
appear to be a local minimum in the rotation velocity due south of the galaxy's center at a distance of i? ~ 1.5i?e (7 kpc). 

Based on the kinematics of dwarf galaxies and X-ray imaging, Binggeli (1999) has argued that the Virgo cluster, and 
particularly its dwarf galaxies, are not yet in dynamical equilibrium. Rather, material appears to be infalling onto M87 
along the principal axis of the Virgo cluster. It is notable that this axis passes though both M87 and M86 — each of 
these supergiant elliptical galaxies is embedded in a swarm of dwarf elliptical galaxies and thereby defines a major Virgo 
subcluster. It is therefore natural to ask if the kinematics of the GCS surrounding M87 is related to that of the more 
distant Virgo galaxies. Such a comparison is especially relevant to the notion that some fraction of the M87 GCS might 
be composed of the accreted and/or infalling remains of "failed dwarfs". 

In order to compare the large-scale kinematics of Virgo directly to that of the M87 GCS, we have carried out the 
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same smoothing procedure described above for the Virgo cluster galaxies. Using the NASA Extragalactic Database, we 
identified a sample of 117 galaxies which are located within 2° (525 kpc) of M87 and have measured radial velocities 
in the range —1000 < cz < 3500 km s~^. The two-dimensional rotation field for this sample is shown in panel (e) of 
Figure [l9|; the corresponding plane of best-fit to the observed velocity residuals is shown in panel (f). Although the limited 
sample size makes firm conclusions impossible, it does appear that the Virgo galaxies show evidence for a gradient in their 
line-of-sight velocity which is roughly in the same direction as that observed for the distant metal-poor GCs surrounding 
M87 — broadly consistent with the infall picture.^ 

Finally, we consider the one- dimensional velocity dispersion profile of the M87 GCS in the larger context of the 
surrounding Virgo cluster. Figure EQ shows the measured and predicted GC aperture dispersion profiles from Figures nT 



15, placed now on an expanded radial scale that also includes data for 246 early-type (mostly dE) Virgo galaxies from 
Girardi et al. (1996). Error bars on the latter measurements refer to 95% confidence limits. The three curves show the 
predicted dap profiles corresponding to the GC subsystems and orbital anisotropics indicated in the lower right corner of 
the figure. Cases of purely radial orbits, /3 = 1, are not shown since all such models are ruled out at high confidence (see 
Figures P^^|l5|).p| Given that the early- type galaxies in Virgo trace the cluster's dark-matter density profile (McLaughlin 
1999a), the striking match seen here between the kinematics of the GCs and the galaxies lends support to our earlier 



assertion (recall eqs. 12 1 |lj] and Fig. |l^) that, in the limit of large galactocentric radius, the M87 GCS also traces the 
Virgo dark matter. Thus, the properties of the globular clusters surrounding M87 are consistent with their being an 
ensemble of test particles — identical in this regard to the Virgo galaxies — which happen to be orbiting in the potential 
well defined by both the galaxy and its parent cluster. 

6. SUMMARY 

We have presented a dynamical analysis of the GC system surrounding M87, the cD galaxy at the dynamical center 
of the Virgo cluster. Our database consists of new wide-field imaging in the Washington C and Ti filters, and radial 
velocities for 278 GCs taken from the literature or obtained during a dedicated spectroscopic survey conducted with the 
Canada- France-Hawaii telescope. This constitutes the largest sample of radial velocities for pure Population II tracers yet 
assembled for any galaxy. A complete description of the database is presented in a companion paper (Hanes et al. 2001). 

Surface density profiles measured from our CTi images reveals the metal-poor GCS to be more spatially extended than 
the metal-rich population, consistent with earlier findings based on HST imaging in the central regions of the galaxy. 
Beyond a radius of i? ~ 1.5i?e (10 kpc), the metal-poor clusters dominate the total GCS. The combined, metal-poor and 
metal-rich samples all show significant rotation of mean amplitude Q,R ~ 170 km s""'^ about axes whose position angles are 
indistinguishable from that of the photometric minor axis, Oq — 65°. The metal-rich GCS shows a roughly fiat rotation 
curve of mean amplitude fii? = 160ltg9° Apart from a possible drift in position angle of the rotation axis at 

R ^ 2Re, the metal-rich clusters appear to be rotating everywhere about the photometric minor axis of the galaxy. The 
rotation field, however, appears non-cylindrical — a two-dimensional map of the rotation field for the metal-rich GCs 
shows a "double-lobed" pattern, with maxima at R ^ 3.5-4i?e (25-30 kpc) along the approximate photometric major axis 
of the galaxy. 

While the mean rotation velocity of the metal-poor GCS, QR — 172^^j!)g km s^^, is indistinguishable from that of its 
metal- rich counterpart, a two-dimensional map of the rotation field suggests that this system may rotate as a solid-body 
or, alternatively, it may exhibit a "shear" in the line-of-sight velocity. This putative shear is similar to that observed 
for Virgo galaxies within two degrees of M87 — broadly consistent with a scenario in which some, or all, of the distant, 
metal-poor GCs are gradually infalling onto M87 along Virgo's principal axis. Inside R ~ 2.5i?e (18 kpc), the approximate 
onset of the galaxy's cD envelope, the metal-poor GCS appears to rotate around the photometric major axis. 

The surface brightness profile of M87 has been combined with ROSAT observations of the surrounding X-ray gas to 
construct a mass model for the Virgo cluster. A comparison of the observed and predicted GC velocity dispersion profiles 
based on this model, and on the measured GC density profiles, suggests that the velocity ellipsoid of the composite GCS 
is almost perfectly isotropic. Dividing the sample on the basis of metallicity reveals the metal-poor clusters to be slightly 
biased to tangential orbits, with /3ci — —0.4, while the metal-rich cluster system shows a radial bias of the roughly the 
same magnitude: /3ci — +0.4. We note that a simple mass model for M87 and Virgo provides a remarkably accurate 
match to the observed GC kinematics, suggesting that the GCs, as a whole, belong to both M87 and the surrounding 
Virgo cluster. More detailed conclusions concerning the formation of M87 and its GCS must await the measurement of 
precise ages for the metal-poor and metal-rich GC systems — the most outstanding observational constraints bearing on 
this issue. 



PC gratefully acknowledges support provided by the Sherman M. Fairchild Foundation during the course of this work. 
DEM acknowledges support from NASA through grant number HF-1097.01-97A awarded by the Space Telescope Science 
Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA under contract 
NAS5-26555. DG acknowledges financial support for this project received from CONICYT through Fondecyt grant 

■^^Note that the observed "shear" is probably not due to the projection of the Virgo space velocity onto the plane of the sky. The expected 
magnitude of this effect is "Uproj = fsSinH, where Vb is the cluster space velocity and Q ~ 0.3° is the angle spanned by the GCs with measured 
velocities. Thus, for an assumed Virgo space velocity of = 2000 km s~^, fproj ~ 10 km s~^, far less than the size of the velocity gradient 
seen in panel (d) of Figure hfl. 

■^^Indeed, the fact that purely radial orbits are so strongly ruled out for the GCS as a whole suggests that any infall component of GCs on 
high-eccentricity orbits probably constitutes only a subsample of the overall system. 



28 



DYNAMICS OF THE M87 GLOBULAR CLUSTER SYSTEM 



8000002, and by the Universidad de Concepcion through research grant No. 99.011.025-1.0. The research of DAH and 
GLHH is supported through grants from the Natural Sciences and Engineering Research Council of Canada. DAH is 
pleased to thank the Directors of the Dominion Astrophysical Observatory and the Anglo-Australian Observatory for 
their hospitality and support during a research sabbatical. This research has made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administration. 

REFERENCES 



Arnaboldi, M., Freeman, K.C., Gerhard, O., Matthis, M., Kudritzki, 
R.P., Mendez, R.H., Capaccioli, M., & Ford, H.C. 1998, ApJ, 507, 
759 

Arp, H. 1968, PASP, 80, 129 

Ashman, K.M., & Zepf, S.E. 1992, ApJ, 384, 50 

Beers, T.C., Flynn, K., & Gebhardt, K. 1990, AJ, 100, 32 

Bahcall, J.N., & Soneira, R.M. 1981, ApJS, 47, 357 

Barnes, J. 1992, ApJ, 393, 484 

Barnes, J., & Efstathiou, G. 1987, ApJ, 319, 575 

Barnes & Hernquist 1996, ApJ, 471, 115 

Bender, R., SagUa, R.P., & Gerhard, O.E. 1994, MNRAS, 269, 785 
BinggeU, B., Sandage, A.R., & Tammann, G.A. 1985, AJ, 90, 1681 
Binggeh, B., Popescu, C.C., & Tammann, G.A. 1993, A&AS, 98, 275 
BinggeU, B. 1999, in The Radio Galaxy Messier 87, Ringberg 

Workshop, ed.H.-J. Roser and K.Meisenheimer (Springer, Berlin), 

9 

Binney, J., &; Tremaine, S. 1987, Galactic Dynamics, (Princeton: 

Princeton Univ. Press), 640 
Carlberg, R.G. 1986, ApJ, 310, 593 
Cohen, J.G., & Ryzhov, A. 1997, ApJ, 486, 230 
Cohen, J.G., Blakeslee, J. P., & Ryzhov, A. 1998, ApJ, 496, 808 
Cohen, J.G. 2000, AJ, 119, 162 
Cote, P. 1999, AJ, 118, 406 

Cote, P., Marzke, R.O., & West, M.J. 1998, ApJ, 501, 554 
Cote, P., Marzke, R.O., West, M.J., & Minniti, D. 2000, ApJ, 533, 
869 

Carter, D., & Dixon, K.L. 1978, AJ, 83, 6 
Davies, R.L., & Birkinshaw, M. 1988, ApJS, 68, 409 
de Vaucouleurs, G., & Nieto, J.-L. 1978, ApJ, 220, 449 
Dubin ski, J- 1998, ApJ, 502 141 



EfstatHiou, G., & Jones, B.J.T. 1979, MMRAS, 186, 133 
FabianI A.C., Nuloon, P.E.J. , &c CaniBaroo, C.R. Nature, 310, 733 



Fall, S.M. 1979, RvMP, 51, 21 
Ferguson, & Binggeli 1994, A&ARv, 6, 67 
Ferrarese, L., et al. 1996, ApJ, 464, 568 
Fillmore, J. A. 1986, AJ, 91, 1096 

Franx, M., lUingworth, G., & de Zeeuw, T. 1991, ApJ, 383, 112 
Frenk, C.S., White, S.D.M., Davis, M., & Efstathiou, G. 1988, ApJ, 
327 507 

Gebhardt, K., Pryor, C, Williams, T.B., & Hesser, J.E. 1995, AJ, 
110, 1699. 

Geisler, D. & Forte, J.C. 1990, AJ, 350, 5 

Geisler, D., Lee, M.G., & Kim, E. 1996, AJ, 111, 1529 

Geisler, D., Lee, M.G., & Kim, E. 2001, in preparation 

Girardi, M., Fadda, D., Giuricin, G., Mardirossian, M., Mezzetti, M., 

& Biviano, A. 1996, ApJ, 457, 61 
Hanes, D.A., Bridges, T.J., Cote, P., Geisler, D., McLaughlin, D.E., 

Hesser, J.E., Harris, G.L.H., & Lee., M.G. 2001, in preparation 

(Paper I) 

Harris, W.E., & Pudritz, R.E. 1994, ApJ, 429, 177 
Harris, W.E., Harris, G.L.H., & McLaughlin, D.E. 1998, AJ, 115, 
1801 

Hernquist, L. 1992, ApJ, 400, 460 
Hernquist, L. 1993, ApJ, 409, 548 

Heyl, J.S., Hernquist, L., & Spergel, D.N. 1996, ApJ, 463, 69 
Huchra, J., & Brodie, J. 1987, AJ93, 779 

Hui, X., Ford, H.C, Freeman, K.C., & Dopita, M.A. 1995, ApJ, 449, 
592 

Jordan, A., Cote, P., West, M.J., & Marzke, R.G. 2002, in 
preparation. 

Kissler-Patig, M., & Gebhardt, K. 1998, AJ, 116, 2237 

Kundu, A., Whitmore, B.C., Sparks, W.B., Macchetto, F.D., Zepf, 

S.E., & Ashman, K.M. 1999, ApJ, 513, 733 
Lee, M.G., Kim, E., & Geisler, D. 1998, AJ, 115, 957 



McLaughlin, D.E. 1995, AJ, 109, 2034 
McLaughlin, D.E. 1999a, ApJ, 512, L9 
McLaughlin, D.E. 1999b, AJ, 117, 2398 

McLaughlin, D.E. 2000, in ESA SP-445, Star Formation from the 
Small to the Large Scale, ed. F. Favata, A. Kaas, and A. Wilson 
(Noordwijk: ESA), p. 77 
McLaughlin, D.E., Harris, W.E., & Hanes, D.A. 1993, ApJ, 409, L45 
McLaughlin, D.E., Harris, W.E., & Hanes, D.A. 1994, ApJ, 422, 486 
Merritt, D., & Tremblay, B. 1993, AJ, 106, 2229 
Merritt, D., & Oh, S.-P. 1997, AJ, 113, 1279 
Merritt, D., Meylan, G., & Mayor, M. 1997, AJ, 114, 1074 
Mihos, J.C, & Hernquist, L. 1994, ApJ, asp437, 47 
Mould, J.R., Oke, J.B., & Nemec, J. 1987, AJ, 93, 53 
Mould, J.R., Oke, J.B., de Zeeuw, P.T., & Nemec, J. 1990. AJ, 99, 
1823 

MuraU, C, & Weinberg, M.D. 1997, MNRAS, 291, 717 
Navarro, J.F., Frenk, C.S., & White, S.D.M. 1997, ApJ, 490, 493 
Neilsen, E.H., Tsvetanov, Z.I., & Ford, H.C. 1999, in The Radio 

Galaxy Messier 87, Ringberg Workshop, ed.H.-J. Roser and 

K.Meisenheimer (Springer, Berlin), 50 

Nulsen, P.E.J. , & Bhringer, H. 1995, MNRAS, 274, 1093 
Ostriker, J. P., Binney, J., & Saha, P. 1989, MNRAS, 241, 849 
Peebles, P.J.E. 1969, ApJ, 155, 393 

Pesce, E., Capuzzo-Dolcetta, R., & Vietri, M. 1992, MNRAS, 254, 
466 

Pierce, M.J., Welch, D.L., McClure, R,D, van den Bergh, S., Racine, 

R., & Stetson, P.B. 1994, Nature, 371, 385 
Quinn, P.J., & Zurek, W.H. 1988, ApJ, 331, 1 
Ratnatunga, K.U., & Bahcall, J.N. 1985, ApJS, 59, 63 



_a 



A. J. 



Kochanek, C. S. 2000, preprint (astro- 



ph/0008062) 

Schlegel, D.J. *^ Finkbeiner, D.P., & Davis, M. 1998, ApJ, 500, 525 
Schombert, J.M. 1988, ApJ, 328, 475 
Schweizer, F. 1980, ApJ, 237, 303 

Seeker, J., Geisler, D., McLaughlin, D.E., & Harris, W.E. 1995, AJ, 
109, 1019 

Strom, S.E., Strom, K.M., Wells, D.C, Forte, J.C, Smith, M.G., & 

Harris, W.E. 1981, ApJ, 245, 416 
Sharpies, R.M., Zepf, S.E., Bridges, T.J., Hanes, D.A., Carter, D., 

Ashman, K.M., & Geisler, D. 1998, AJ, 115, 2337 
Tonry, J.L. 1987, in The Structure and Dynamics of Elliptical 

Galaxies, lAU Symposium No. 117, ed. T. de Zeeuw, (Dordrecht, 

Reidel), 89 

van der Marel, R.P. 1991, MNRAS, 253, 710 
van der Marel, R.P. 1994, MNRAS, 270, 271 

Wahba, G. 1990, Spline Models for Observational Data (SIAM, 
Philadelphia) 

Warren, M.S., Quinn, P.J., Salmon, J.K., & Zurek, W.H. 1992, ApJ, 
399, 405 

Weil, M.L., & Hernquist, L. 1994, ApJ, 431, L79 

Weil, M.L., & Hernquist, L. 1996, ApJ, 460, 101 

Weil, M.L., Bland-Hawthorn, J., & Malin, D.F. 1997, ApJ, 490, 664 

West, M.J. 1993, MNRAS, 265, 755 

West, M.J., Cote, P., Jones, C, Forman, W., & Marzke, R.O. 1995, 

ApJ, 453, L77 
West, M.J., & Blakeslee, J. P. 2000, ApJ, 543, L27 
White, R.E. 1987, MNRAS, 227, 185 

White, S.D.M. 1982, in Morphology and Dynamics of Galaxies, 
Twelfth Advanced Course of the Swiss Society of Astronomy 
and Astrophysics, ed. L. Martinet & M. Mayor (Geneva: Geneva 
Observatory), 291 
Zepf, S.E., & Ashman, K.M. 1993, MNRAS, 264, 611 
Zepf, S.E., Beasley, M.A., Bridges, T.J., Hanes, D.A., Sharpies, R.M., 
Ashman, K.M., & Geisler, D. 2000, AJ, 120, 2928 



This figure "m87paper2-figl8.jpg" is available in "jpg" format from: 



http://arXiv.0rg/ps/astro-ph/0 1 06005 V 1 



